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1 Why do we need lattice QCD? 



The theme of this year's TASI is CP violation. CP violation offers a possible window 
on physics beyond the Standard Model: Will the single phase in the CKM matrix 
be sufficient to explain all the CP violating amplitudes that we hope to measure 
once B-factories, </>— factories and the next generation of e'/e experiments are up 
and running? To answer this we need to relate the phases in the CKM matrix, 
which appear in the underlying quark amplitudes, to the measurable phases in 
hadronic decay amplitudes. For some B-meson decays, e.g. B — > tpKs, we can avoid 
uncertainties due to hadronic structure by taking appropriate ratios. For most 
quantities, however, the relation between CKM and measurable phases is obscured 
by non-perturbative physics. To find the relation, we need to know certain hadronic 
matrix elements of fermion bilinears and quadrilinears. One of the most practical 
uses of lattice QCD is the calculation of such matrix elements. 

Figure |I| s hows a cartoon of one of the best studied examples, CP-violation in 
the — K'^ mixing amplitude. Most of the attendees of TASI were not born 
when this amplitude, parameterized by e, was first measured. We understand the 
central "box" in the diagram: the large mass of the top quark allows us to treat it 
as an effective four-fermion operator, with a coefficient proportional to Im[Vj^Vj^], 
multiplied by a calculable QCD anomalous dimension factor. What we do not know 
how to calculate analytically is the effect of low momentum gluons and quarks 
{\p\ < 2 GeV, say). Such gluons confine the q-q pairs into kaons. They interact 
with a large coupling constant, as{p), and thus the interactions cannot be described 
using perturbation theory. Non-perturbative calculations are needed. The only 
method presently available which starts from first principles and makes no further 
assumptions is to simulate lattice QCD (LQCD) numerically. In the example of e, 
what the lattice must provide is the matrix element 
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Here the decay constant is normalized so that /^r = 93 MeV. 
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Figure 1: The K — K mixing amplitude. 



I will discuss the lattice calculation of Bk, which is just a parameterization of 
the matrix element, at some length below. What I want to point out here is that 
one needs to know in order to use the experimental result for e to extract a 
value for the Im(V^^). A similar matrix element is needed to extract information 
from the measured B — B mixing amplitude. Other matrix elements are needed to 
extract information from semileptonic S-meson decays. 

Thus to probe the electroweak theory one must be able to do non-perturbative 
calculations of matrix elements. These lectures concern only lattice calculations 
of these matrix elements. It is important to realize, however, that there are other 
methods available: the large (number of colors) approach, QCD sum rules, 
chiral quark model, etc. These all have the status of (more or less sophisticated) 
models: they make approximations which are hard to improve upon and whose 
effect is not always easy to gauge. They often have the advantage, however, of 
being applicable to a wider range of quantities than can be studied on the lattice.^ 
What has happened in the last few years is that lattice calculations have come of 
age, in the sense that for some quantities all errors are understood and are small. 
Future progress will, I think, see the lattice give reliable answers for increasingly 
many quantities. This will allow us to test and improve the approximate methods, 
which can then be applied with more confidence to quantities for which lattice 
calculations are difficult. 

Lattice calculations involve, at their core, numerical simulations, but in order to 
make the various extrapolations that are needed, considerable guidance from ana- 
lytic results is required. Thus the lattice phenomenologist needs, in her or his tool 
kit, not only a PSC (personal super-computer), but also expertise in lattice per- 
turbation theory (to match lattice onto continuum operators), chiral perturbation 
theory (to extrapolate to physical quark masses), non-relativistic QCD or its close 
variants (to study heavy quarks on the lattice), and the technology of "improved 
actions" (how to reduce errors due to finite lattice spacing). In the following I hope 
to give a flavor of how these tools are used, and to indicate the present status of 
their application. 

2 Basics of Euclidean Lattice Field Theory 

I will begin with a review of the basics. My discussion will be both sketchy and 
patchy. Those wishing more detail or rigor will find both in two recent texts on 
lattice field theory [|1], 0. Less detailed but still useful is Creutz's monograph 0. 
Many results I use are standard, and if I do not give a reference, it can be found in 
one or more of these books. 

The steps that are taken to get to a numerical simulation are these: 



*I will discuss the limitations of lattice calculations in the following. 
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1. Use the Euclidean space functional integral formulation of field theory 

Z = y [dA][dq][dq]e- , (2) 

where A, q and q are respectively the gauge, quark and antiquark fields. 

2. Discretize the theory on a hypercubic lattice, with spacing a. In present 
simulations this lies in the range 0.05 — 0.2 fm. 

3. Work in finite volume: L points in the three spatial directions, T points in 
the time direction. The largest lattices in use today are roughly 32^ x 64 in 
size, which, for a = O.lfm, corresponds to (3.2 fm)^ x 6.4 fm. 

4. Do the functional integral (now a multidimensional path integral) using nu- 
merical Monte Carlo methods. 

5. Attempt to extrapolate L ^ oo and a ^ 0. 

Actually this is something of an idealization. Additional steps are required in 
most simulations. 

3.5 Make the "quenched" approximation, in which internal quark loops are left 
out, keeping only valence quarks. 

5.5 Extrapolate from quark masses ruq ~ ms/2 down to the physical up and down 
quark masses (m„ + md)/2 ~ ms/2h. 

Many of these steps will be discussed in more detail below. 

It is important to realize that, with the exception of "quenching", the approx- 
imations that are made can be systematically improved. This improvement can 
occur not only because of increases in computer power (which allows one to study 
larger lattices, for example), but also due to analytical advances (e.g. "improving 
the action" so that one can work with larger lattice spacings without increasing the 
errors due to discretization). 

An important limitation of LQCD is that the calculations are carried out in 
Euclidean space. Why is this? Because the integrand in a Minkowski path integral, 
exp(iS'M), is complex. In contrast, the Euclidean integrand is, in most cases, real 
and positive. Thus in the Minkowski integral there are large cancellations between 
different regions of configuration space, and these make it hard to simulate all but 
very small systems. This is an algorithmic, not a fundamental, issue. It is possible 
that it will be resolved in the future, though there are no signs of this at present. It 
turns out that even in Euclidean space, QCD at finite chemical potential (i.e. finite 
baryon number density) has a complex action, and thus is very difficult to simulate. 
Similarly, chiral theories have complex Euclidean actions. In fact, the "fermion 
doubling problem" (to be discussed below) makes it difficult to even formulate such 
theories on the lattice. 
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2.1 Extracting information from Euclidean Correlators 

In this subsection I want to explain why LQCD is well suited to the calculation of 
matrix elements involving the vacuum or single particle states, while those involving 
multiple particles are much more difficult. What can be most successfully calculated 
are 

• The energies of states, e.g. (7r(p)|if|7r(p)), which, in the continuum limit 
should become E = ^jP + m^. In particular, for p = 0, one directly calculates 
particle masses. 

• Decay constants, e.g. (OIA^Itt) = i\/2p^f^^ where is the axial current. 

• Single particle matrix elements {N{p)\0\M{q)), where and M are single 
particle states, and (9 is a fermion bilinear, or quadrilinear, or perhaps a 
gluonic operator. 

I will explain why these quantities are easily calculable by showing how they are 
calculated. Consider the Euclidean two-point function 

C(r) = J[dfi]exp{-SE)O*iT)Oi0) (3) 
= (O*(r)O(0)) (4) 

where is the measure for the quark and gluon fields, Z = /[rf/i] exp(— S'^) is 
the partition function, and 0{t) is an function of the fields residing at Euclidean 
time r. To study a pion at zero three-momentum, for example, we might choose 
O = J2x^{t, x)75d(r, x). Recall that the fermion fields in the functional integral are 
Grassman variables. 

The functional integral is constructed to give a time-ordered expectation value 

C{r) = {0\T[d\T)dm\Q) (5) 
= (O|e^"ate-^"0|O) (6) 
= (OlOV^^aiO) , (7) 

where H is the Hamiltonian operator, normalized so that its ground state has zero 
energy, H\0) = 0, and 0{t) is the Heisenberg operator corresponding to C(r) 
(in simple cases obtained by substituting for each field in O the corresponding 
operator). In the second line I assume r > 0, so that r-ordering, T, has no effect, 
and I have used the Euclidean version of time translation to shift the operator back 
to r = 0. By convention O{t = 0) = O. 

Before proceeding, it is worthwhile understanding the conditions under which 
expectation values in Euclidean functional integrals (e.g. C(r)) can be written as 
time-ordered products, as in Eq. |[ Textbooks typically begin with the Minkowski 
space time-ordered product, which is then analytically continued to Euclidean space 
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(in the form of Eq. |^), and then written as a functional integral. What we want to 
do is run the argument the other way around. This question was studied long ago by 
Osterwalder and Schrader, who found the following If the action Se is Euclidean 
invariant, and expectation values satisfy a property called "reflection positivity", 
plus some other more technical conditions, then there exists a Hilbert space with 
positive norm, a Hamiltonian H acting on this space which is hermitian, and whose 
spectrum is bounded from below with the lowest state having zero energy, and fleld 
operators, such that Eqs. hold.|] There is then no obstruction to analytically 
continuing these correlators to complex time by an inverse Wick-rotation, r — > r' = 
e^'^T. For = 7r/2, what results are Minkowski time-ordered products 

CiT = tt) = {0\T[d\t)dm\0) , (8) 

where t is Minkowski time. A consequence of the initial Euclidean invariance is that 
there are unitary operators acting on the Hilbert space which implement Poincare 
transformations. 

Once we have these time-ordered products we can, in principle, use the LSZ 
reduction formalism to extract the S-matrix from the residues of their poles. In other 
words, Euclidean functional integrals with suitable properties provide a deflnition 
of the fleld theory, as long as we can carry out the analytic continuations. It is 
important to realize that not all Euclidean functional integrals satisfy the necessary 
conditions. In particular, QCD in the quenched approximation can be written as 
a particular functional integral, displayed below, but does not satisfy reflection 
positivity, and does not correspond to a well-behaved Minkowski fleld theory. 

With this background, let us return to the correlator C(r). Inserting a complete 
set of energy eigenstates (single and multiparticle states), we flnd 

C(r) = El(n|O|0)p|— -. (9) 

Here I am assuming a flnite volume V , and using relativistically normalized states, 
m = ^EV5,^,J,^,^6,^,^ 2Ei2nf5'ip- q) , (10) 

where = + m^. We could now follow the LSZ procedure: Fourier transform 
to Euclidean energy, rotate to Minkowski space, and look for poles. We would flnd 
a series of poles corresponding to the stable states which couple to the operator 
O, and various cuts for multiparticle intermediate states. If the lightest state is a 
stable particle, however, then there is no need to rotate to Minkowski space. For 
example, if the operator creates a tt'^ at rest, as in the example above, then we can 
read off simply by looking at the exponential fall-off at large Euclidean times 

C{r) |(7r+(p = 0)|a|0)pexp(-m.r) . (11) 
tPor a nice discussion, including the definition of reflection positivity, see Ref. 0. 
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Figure 2: "Kaon" correlator with staggered fermions. 

Furthermore, the coefficient of the exponential gives us a vacuum to pion creation 
amplitude, which, if O = u'jQ'j^d, is proportional to the decay constant For 
this procedure to work it is important that there be a gap between the lightest and 
next-to-lightest states. In the present example, the latter consists of three pions at 
rest, with E = 3m^, so there is a gap for finite pion mass. 

Figure ^ shows an example of a two point function computed numerically (from 
a calculation done in collaboration with Greg Kilcup and Rajan Gupta). It is for a 
particle with the quantum numbers of the kaon, and with a mass similar to that of 
the kaon, but in which m^"* = m^'** ^ ^m^. The lattice spacing is about a = 1/12 
fm, and the lattice size is 32^ x 48. The graph shows that by r ~ 20, the data are 
represented almost perfectly by a single exponential. The solid line shows a fit to 
such a form in the range r = 20 — 37. Outside this range, the dashed line shows 
the extension of the fit function. Thus you can see the deviation of the data from a 
pure exponential at short times. The curvature for r > 40 is due to the boundary 
conditions. 

The observant reader will notice an inconsistency between the data and the 
expected form, Eq. ^. C(r) is a sum of exponentials with positive coefficients, and 
must approach its asymptotic form from above. This condition does not apply to 
the results of Fig. |^, however, because we use different operators at times and 
t. This allows the coefficients to have either sign, and the approach to asymptotia 
need not be monotonic. One of the technical aspects of lattice calculations which 
I will not go into, is the use of improved operators (sometimes called "sources"), 
designed so as to couple strongly to the lightest state but much more weakly to 
higher states. With such operators the correlators quickly become dominated by a 
single exponential. This allows one to use lattices which are shorter in the Euclidean 
time direction, and usually reduces the signal to noise ratio. 
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Figure 3: Example of an effective mass plot. The fitting range is shown by dotted 
lines, the fit value by the solid line. 



It is common for lattice practitioners to present their results in terms of the 
logarithmic derivative —dln[C{T)]/dT, which, for r — > oo, becomes the energy of 
the lightest state. The lattice version of this is arues = ln[C(r)/C(r + l)], where a is 
the lattice spacing. An example of such a plot is shown in Fig. ^ for a calculation 
at similar parameters to that in Fig. ^ (more precisely, a 30 x 32^ x 40 lattice at 
/3 = 6.17, K = 0.1532), taken from Ref. 0. (The behavior for r > 20 is due to a 
negative parity j = | particle propagating backwards in time.) Note that, for the 
operators used to obtain these results, the coefficient of the non-leading exponential 
turns out to be positive, as shown by the fact that rrics approaches its asymptote 
from above. 

To first approximation the mass can just be read off from the graph. It is not 
possible, to use the plot to give an idea of how well the data is represented by a 
single exponential. This is because all the points are correlated, fiuctuating up and 
down nearly in unison. One needs the full correlation matrix, and not just the 
"diagonal" errors which are displayed, in order to test the goodness of fit. It is true, 
however, that the fiuctuations increase with r, so one is always balancing the need 
to go to longer times, so as to be sure one has a pure exponential, with the desire 
to work in the region where the errors are smaller. This is a general feature of the 
analysis of lattice "data" . 

Single particle matrix elements can also be calculated directly in Euclidean 
space, e.g. 

C(r2, n) = {U{p, t^)OkM ri)K{0)) . (12) 
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Here, the operators K ~ ^755 and 11 ~ U'~^^d create a and destroy a 7r~, re- 
spectively. (I am using the loose description in which the I associate the functions 
appearing as arguments in the functional integral with the corresponding operators 
which appear in the operator representation.) The operator in the "middle" (as- 
suming r2 > Ti > 0), might be, for example, the current 57^(1 +75)^. Using the 
generalization of Eq. ^ one finds 

C(r2,ri) = (0|n(p)e-^(^^-^^)O^^(g-)e-^^i]r|0) (13) 
= Y.m{p)\n) {n\OKM\n')^^{n'\m . (14) 

If T2 — Ti and Ti are both large enough, then we need keep only the lightest state in 
the two sums. Thus {n\ = (7r°(p)|, \n') = \K^{k), where k = p — q. The correlator 
becomes 

C{T2,n) oc {0\U{p)\7r\p))e-^-^^^-^^H7r\p)\6^^\K\k^^^^ . (15) 

The energies and the creation and annihilation amplitudes can all be obtained from 
two point correlators, and divided out. Thus, from the three-point function one 
can directly extract the transition matrix element {tt^ {p)\0 Kn\K^ {p}) ■ If C ~ ^7^ti, 
then this is the vector form factor which governs the semileptonic decay spectrum 
in K*^ -K^e^v. By changing the quantum numbers of the operators one can 
calculate other form factors of interest, e.g. B K*'y. Using a four-fermion 
operator for Ok^ one can calculate K — K and B — B mixing amplitudes. Most of 
the phenomenologically interesting results from LQCD are for matrix elements of 
this type. 

This exhausts the types of quantity that are simple to calculate in Euclidean 
space. It is much more difficult, unfortunately, to calculate amplitudes involving 
two or more particles in either initial or final state. Some of the quantities we would 
like to calculate are 

• ^(tttt — > vttt), 

• A{K — > vrvr), and 

• A{B^i)K,). 

The difficulty arise from the fact that, in Minkowski space, these amplitudes are 
complex. This follows from unitarity, as there are on-shell real intermediate states. 
(The only exception is the pion scattering amplitude at threshold, for which there 
is no phase space for intermediate states.) But on-shell intermediate states are only 
possible in Minkowski space; there are none if the external states have Euclidean 
momenta. Starting from Euclidean momenta, the imaginary parts are generated 
hy the analytic continuation to Minkowski space. A simple example is the fuction 
ln(4m^ — (pi +^2)^), real in Euclidean space, but imaginary upon continuation to 
physical momenta satisfying (pi +^2)^ < — 4m^. 
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One way of seeing why there is no simple way of doing the calculation directly in 
Euclidean space is the following. Consider the K ^ ttti amplitude. This is obtained 
by creating the kaon, acting with the weak Hamiltonian to turn it into a state 
with the quantum numbers of two pions, and then destroying the two pions. The 
physical amplitude involves the two pions having a non-zero relative momentum. In 
a Euclidean correlator, however, one does not obtain this contribution by making 
a large Euclidean time separation between the weak Hamiltonian and the pion 
operators. Instead, what dominates is the transition from a kaon to two pions at 
rest, the latter being the lowest energy state. Even if one uses two pion operators 
having relative momentum, they will have, due to interactions, a coupling to the 
lowest energy state. Thus what dominates the Euclidean correlator is an off-shell 
transition amplitude {K\T-Cw\t^t^{p = 0)), where p is the relative momentum, and 
Tiw the weak Hamiltonian. This is not the quantity of interest. Another way of 
stating the problem is that one does not create the in and out states directly in 
Euclidean space. See Ref. for a clear explanation of this point. 

Thus to get the correct amplitude, both in magnitude and in phase, one has to 
analytically continue. In most cases this will only be possible if one has a model of 
the momentum dependence of the amplitude. For example, for K tttt decays, one 
use chiral perturbation theory, which, at leading order, relates the decay amplitudes 
to calculable single particle matrix elements. Using such a method, however, one 
gives up on the possibility of a first-principles calculation, the errors in which can be 
systematically reduced. There will be an irreducible error due to the uncertainties 
in the model used. For the example of K tttt decays, the uncertainties are due 
to higher order terms in the chiral expansion. 

Another approach is possible in the case of scattering. Liischer has shown, in 
a beautiful series of papers [0, how to extract scattering amplitudes indirectly, 
using the finite volume dependence of the energies of two particle states. Here too, 
however, one must make approximations to use the results in practice. An infinite 
tower of partial waves contribute to two particle energy shifts, and one must assume 
that only a finite number are important. 
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3 Discretizing QCD 



To calculate amplitudes numerically, we must discretize QCD. This must be done 
in such a way that gauge invariance is maintained, since this invariance is required 
to guarantee the unitarity of the S-matrix. How to do this was worked out long ago 
by Wilson. 

3.1 Continuum QCD, a brief overview 

The continuum action is given by 

Se^- E [tP + + i / Tr(F^.F^.) , (16) 

q=u,d,s,c,b,...^ ■'^ 

where the integrals run over Euclidean space. The covariant derivative is 

D^^d^- igA^ , (17) 

in which the gauge fields are collected into a matrix — A^T"', with the 
generators of the SU{3) Lie Algebra 

[T«, T^] = if^bcj^c ^ tr(T«T'') = ISab ■ (18) 

The quark fields are color triplets, with an implicit color index. Finally, the gauge 
field strength is 

F^, = F;,T- = ^p^, D,] = d^A, - d,A^ - ig[A^, A,] . (19) 

A local SU{3) gauge transformations is described by a space-time dependent 
element V{x) e SU{3) (V-^ = V\ (\et{V) = 1): 

q{x) V{x)q{x) q{x) -^q{x)V-\x) (20) 
A^{x) ^ V{x)A^{x)V-\x) + iV{x)d^V-\x) (21) 

F^,{x) ^ V{x)F^,{x)V-\x) (22) 
[D,q]{x)^V{x)[D,q]{x) . (23) 

Given the last two lines, it is simple to see that Se is invariant. 
It is useful to introduce the path-ordered integrals 

L(x,y) ^Pexp{ig [ dz^Ai,{z)} , (24) 
■'y 

which are to be thought of as going from y to x. The ordering is such that, for 
example, A^i^x) is always to the left of An{y). Gauge transformation properties 
depend only on the end points of L, and not on the path of integration 

L{x,y)^V{x)L{x,y)V-\y) . (25) 
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They thus transport the gauge rotation from one point to another, such that the 
quantity q{x)L{x, y)q{y) is gauge invariant. Another gauge invariant quantity is the 
trace of the path-ordered integral around any closed path 

Tr[L(x, x)] — ^Tr[V{x)L{x, x)V-^{x)] = Tr[L{x, x)] (26) 

These objects are called Wilson loops. 

With these quantities in hand, we can now construct a gauge invariant lattice 
version of QCD. One cannot simply place the quark and gauge fields on the sites 
of the lattice and discretize the derivatives appearing in Se- Instead, the gauge 
fields, which transmit information about gauge transformations from one position 
to another, will live on the "links" or "bonds" connecting the sites. I will choose 
the lattice to be hypercubical, since this is the form most easily studied numerically, 
and nearly all work has been done with it. The notation for the sites and links on 
the lattice is shown in Fig. ^ 

Discretizing fermions presents its own set of problems, not directly related with 
gauge invariance. Thus I consider first the gauge part of the action. This will be 
constructed from elements of SU{3), Un,^, associated with the link from site n to 
site n + fi, and corresponding to the continuum line integral along the link: 



Un,^ ~ L(an, an + ajl) 

rn 

= Pexp ig / dZf,Af,{z) 

Jn+afi 

= 1 — igaA^{n+^fi) + 0{a'^) . 



(27) 



The "~" in the first line means "corresponds to" . The vagueness here is deliberate — 
once we put the theory on the lattice there are no gauge fields A^: they are replaced 
by the t/'s. The expansion in the last line is useful, however, for thinking about 
what the t/'s mean, and also for taking the classical continuum limit. 
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The gauge transformation properties of the t/'s are taken to be the same as 
those of the corresponding L's 

Un,^_L VnUn,f_iVn+^ , (28) 

where the Vn G SU (3) are gauge transformation matrices which hve on sites. I have 
adopted the abbreviated notation in which n + /i means n + afi. In correspondence 
with the continuum result L{x, y) = L{y, xY, we associate [/^^ to the hnk from 
n + fi to n. Note that we can multiply the t/'s along any closed loop and take the 
trace, and obtain an object which is invariant under gauge transformations, since 
VnV^ = 1. These are the lattice versions of Wilson loops. 

We can construct a lattice version of the pure gauge action using the smallest 
Wilson loop, that around an elementary square or "plaquette" 

Pi = Ur^.,,Ur,^,,.uUM,u ■ (29) 
The geometry is illustrated here. 



It is reasonable that such a loop is related to F^jy, because the field strength is the 
curvature associated with the connection A^. In any case, using the correspon- 
dence given above for the ?7's, and after some algebra, one finds that the classical 
continuum limit of the plaquette is 

Pi = 1 - tga^F^, - 'y^^FI + la^'G^, + la^H^, + O(a^), (30) 

where H^, and G^, and are hermitianQ. Thus one can use the iiv plaquette as a 
discretized version of the corresponding component of the field strength, F^j^. If we 
take the trace, so as to get a gauge invariant quantity, we find 

Re TrP^, = N,- iaHr{Fl) + 0(0*5) ^ (g^) 
where iVc = 3 is the number of colors. We then have 

/ d\ J2 iTrF^-F^. ~ Jl^i^c- ReTrD) . (32) 

The factor of 2 arises because of the mismatch between the number of plaquettes 
per site, 6, and the number of terms in the sum J^/iu^ 12- 

■^Exercise: derive this result. Everyone should do it once! 
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Now, it is standard (though unfortunate) to replace the couphng constant by 
^ = so that 

Sg = -PY, + irrelevant constant . (33) 

□ 

This is called the "Wilson (gauge) action". It is important to realize that there 
is nothing special about using the smallest loop to define the action. Any loop, 
e.g. a 1 X 2 rectangle, contains a term in its expansion proportional to a linear 
combination of components of (F^,y)^. By taking an appropriate combination of 
loops we can obtain the continuum action as a ^ 0. The advantage of a small loop 
is that corrections proportional to powers of the lattice spacing are typically smaller 
than with a larger loop. 

At this juncture it is appropriate to make a brief comment on "improving" the 
action, i.e. reducing the errors due to discretization. In Eq. ^ there are correc- 
tions of O(a^) compared to the term that we want. Symanzik has shown how to 
systematically reduce these corrections from O(a^) to 0{g^a?) (by one loop calcula- 
tions), and then to 0{g'^a^) (by two loop calculations, etc. p. More ambitious is the 
program (based on using the renormalization group) to construct an almost perfect 
action, i.e. one in which all terms of 0(a^") are almost absent [0. The subject has 
lain dormant for almost a decade, but is now receiving considerable attention, in 
part because progress has been made at reducing the errors in the fermion action 
(which are usually of 0(a), and thus larger than those in the gauge action). To 
date, however, numerical simulations leading to phenomenological results have only 
been carried out using the Wilson gauge action. 

To complete the definition of the theory I need to specify the measure. Each 
link variable is integrated with the Haar measure over the group manifold. This 
measure satisfies {V and W are arbitrary group elements) 

j dUF{U) = J dUFiUV) = J dUF{WU) . (34) 

Given this, it is simple to see that the functional integral 

^gauge = / n dUn,^, ^Xp ( E 1^^^ TrD ) (35) 
links \ □ " / 

is gauge invariant. What has been accomplished here is a non-perturbative, gauge 
invariant regularization of pure gauge theories. What has been sacrificed is full Eu- 
clidean invariance: rotations and translations. The hope is that, as one approaches 
the continuum limit, these symmetries are restored. 

Pure SU{2i) gauge theory is not QCD. But it is still of interest as a non-trivial 
field theory, sharing some properties in common with QCD. In particular, its spec- 
trum should consist of massive glueballs, in which the gluons are confined by their 
self-interactions. Considerable progress has been made in numerically simulating 
this theory. I will give more details below of the methods used. For now, let me 
note that one calculates the glueball spectrum by looking at two point functions 
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I 1 

Figure 5: Spectrum of pure gauge SU(3) theory. 

in which the operators are Wilson loops of various shapes and sizes. By forming 
appropriate linear combinations one can project onto the various representations of 
the lattice rotation group, which can then be associated with certain spin parities 
in the continuum limit. The present status of the spectrum is shown in Fig. ^ 
(Ref. [0). The masses are given in units of the square-root of the string tension, 
a quantity I discuss below. For the moment, just consider the units as arbitrary, so 
that all we can extract is the ratio of the masses of different glueballs. Clearly, we 
have reasonable control over the spectrum of this non-trivial theory, with the scalar 
glueball being the lightest, followed by the tensor and pseudoscalar. Furthermore, 
there is evidence that Euclidean symmetry is being restored. For example, the five 
spin components of the tensor glueball lie in two different representations of the 
lattice cubic rotation group, yet all have the same mass within errors. 

I have not explained yet how one takes the continuum limit. When one does a 
simulation, one picks the value of /3 = 6/(7^, not the lattice spacing. The output of 
the calculation are a set of masses in lattice units: amgiuc- One obtains a value for 
a by comparing these to a physically measured scale. In this case, there is no such 
scale, since the real world is not pure gauge theory. If it were, we would use the 
physical value of mgiuc, for a particular glueball, to fix a. All other masses are then 
predictions. If we choose a different value for (3 (i.e. for g'^), then we would obtain 
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a different value for a. To approach the continuum hmit we adjust (3 so that a —> 0. 
As we do this, the lattice mass vanishes, anigi^e 0. Thus correlators fall off more 
and more slowly, C{t) oc exp(— amgiue^)- It is conventional in statistical mechanics 
to describe this in terms of the divergence of the correlation length in lattice units: 
C{t) oc exp(— t/^), = l/(amgiue) ^ oo. A place at which ^ diverges is called a 
critical point, and it is at such points that one can take a continuum limit. 

In the case of pure gauge theories, or QCD with sufficiently few flavors that 
it is asymptotically free, one knows that the critical point is at = 0, j3 = oo. 
This is because the functional dependence of on a is calculable for small g^ using 
perturbation theory, yielding the familiar result that (7^ oc l/ln(l/a). Inverting 
this, we flnd that a = (1/A) exp(— 1/2/50(7^) (/?o the leading order coefficient in the 
beta-function), so that a — ^ for (7^ ^ 0. What we do not know a priori is the 
proportionality constant 1/A. This we must determine numerically. Given A, we 
can then predict the variation of a with g'^, for small enough g^. This works well for 



/3 > 6, as long as one chooses an appropriate scheme for the coupling constant [11 



3.2 Discretizing Fermions 

Discretizing the fermionic action, 

Sp = -^ipE + m)^ , (36) 

seems straightforward at first sight. The subscript on the derivative is a reminder, 
immediately to be dropped, that we are in Euclidean space, with hermitian gamma 
matrices satisfying = 25^^ . To discretize Sp, we place quark and antiquark 

fields on sites, q{x) — > and perform gauge transformations like so: 

qn Vnqn , In ^ InVl ■ (37) 

To obtain a discrete derivative we must separate the q and q fields, and we make 
this gauge invariant using the link variables. For example, choosing a symmetric 
derivative on the lattice, 

qD^q{x) — > I g„ [f^n,/.gn+|. - f^I-^.^gn-^.] ■ (38) 

Here and in the following I have set the lattice spacing a to unity. It is left as a 
simple exercise to show that when one expands out the f/'s in terms of ^'s one 
indeed finds the covariant derivative. With this choice of the action is 

- SN{q, q,U) = J2 l^nlf^ [f^n,M?n+M " Ul-^^^q^-f,] + J2 "^^n?n • (39) 
nfi n 

For reasons about to be described, this is called the "naive" fermion action. The 
full partition function of QCD is then 



Zqcd = / n ^^".A* n [dqdq] exp 



links q=u,d,s, 



-SgiU)- SN{q,q,U) 



q=u,d,s,. 



(40) 
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The quark measure is gauge invariant because a special unitary rotation has unit 
jacobian. Thus Zqqb is gauge invariant. 

There are many choices available when discretizing derivatives aside from that 
of Eq. the forward derivative (9^g g„+^ — the backwards derivative 
{d^q qn — qn-^l)^ etc. The symmetric derivative is, however, the most local choice 
which preserves the anti-hermitian nature of the continuum operator p. 

The harmless looking action of Eq. 3£ gives rise to an infamous problem: in 
d dimensions, it represents 2"^ degenerate Dirac fermions, rather than one. This 
sixteenfold replication in 4 dimensions is referred to as the "doubling (!) problem". 
Doubling has nothing to do with gauge fields, and so to discuss it I consider a free 
lattice fermion. As discussed above, the spectrum of the states can be determined by 
looking at the fall-off of the Euclidean two point function. More useful for analytic 
studies is to transform to momentum space, and look for poles. These occur at 
complex Euclidean momenta: k"^ = — m^, where m is the mass of the state. 

To evaluate the two point function, i.e. the propagator, recall the integration 
formulae for Euclidean fermions (which are Grassman variables): 

Z = y[dg][rfg] e^(^+'")'' = det(|) + m) , (41) 
G{x,y) = -Z-^ j[dq][dqY^^^^^'^q{x)q{y) = . (42) 

To diagonalize p we go to momentum space. Due to the periodicity of the lattice 
one is restricted to the first Brillouin zone 

Qn = £ (0e^'"g(fc) , qn = #e-^-g(fc) , (43) 
In terms of these fields, and using = sinfc^, 

-Sn= fqik)iiJ2 s^l^ + m)q{k) . (44) 
Thus, the momentum space propagator is given by 

G{k) = = / , . 45 

ip + m + 

It is useful to reinstate factors of a, so that m = aruphys and k = akphys- In the 
continuum limit, for fixed physical quark mass, m ^ 0. There is thus a pole near 
k = 0, and we can expand = ak^p^ysil + 0(a^)), yielding 

aC{k) = . (46) 

phys ~^ ^^phys 

This has a pole at k^f^y^ = —rn^phy^, representing the fermion that we expected to 
find. 
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Now we come to doubling. The lattice momentum function vanishes for 
fc^ = vr as well as = 0. In the neighborhood of the momentum (tt, 0,0,0), if we 
define new variables hj k[ = n — ki, k[ = ki, i = 2 — 4, then 

To bring the propagator into the standard continuum form, I have introduced new 
gamma- matrices, 7^ = —71, 7- = 7^, z = 2 — 4, unitarily equivalent to the standard 



set. Equation ^ shows that there is a second pole, at fc'^ = — m^, which also 
represents a continuum fermion. This is our first "doubler" . 

The saga continues in an obvious way. vanishes if each of the four components 
of k^ equals or tt. There is a pole near each of these 16 possible positions. Our 
single lattice fermion turns out to represent 16 degenerate states. 

It is not, in fact, the replication of fermions which is the hard part of the problem, 
but rather the way in which the chiralities of the states work out. If m = 0, then 
we can introduce a chiral projection into the action 7^ —>■ 7^(1 + 75)/2, which in 
the continuum restricts one to left-handed (LH) fields. On the lattice, the pole 
near = is LH. The second pole I uncovered, however, actually represents a RH 
field. This is plausible, because 75 = —75 and so (1 + 75) = (1 — 75). It can be 
demonstrated by considering the coupling to external currents. For each of the 
components of k that is near to n the chirality fiips, so that one ends up with eight 
LH and eight RH fermions. It is important to note that, when one introduces gauge 
fields, all the fermions are necessarily coupled in the same way. Thus one always 
obtains a "vector" representation of fermions, i.e. one in which LH and RH fields 
lie in the same representation of the gauge group. 

How general is this result? Karsten and Smit showed that, in infinite volume, 
any local, antihermitian discretization which maintains translation invariance will 
give rise to LH and RH fermions in pairs There need be only one such pair: 
Wilczek has given an example, using an action which breaks Euclidean rotation 
T^ . What are the consequences of Karsten and Smit's result? 



mvariance 



Lattice regularization automatically takes care of the fact that theories with 
anomalous representations of fermions (e.g. SU (Nc) with a single left-handed 
fermion) cannot be defined. 

It does too good a job, however. One cannot discretize a chiral theory, i.e. 
one having an anomaly free, yet chiral, fermion representation. In particular, 
one cannot discretize the electroweak theory. 

Indeed, one cannot even discretize QCD with (nj) massless quarks, in the 
following sense. Such a theory should have an SUijif) l x SU {nf)^ chiral sym- 
metry, under which the LH and RH quarks rotate with independent phases. 
But the lattice fermions are all begotten of the same lattice field, and so cannot 
be rotated independently. 
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In summary, then, the lattice theory lacks chiral symmetries. 

Can this be fixed up? Much effort has been devoted to this question. There are 
various escapes, each with their own problems. Most notable are these. 

• One can explicitly break chiral symmetry right from the start, and aim to 
recover it only in the continuum limit. This is, after all, what one does with 
the rotations and translations. For fermions in vector representations, this is 
the approach originally taken by Wilson, which I discuss in more detail below. 
For chiral theories, this is the approach advocated by the Rome group, and 



involves breaking the gauge symmetry at finite lattice spacing [14|. 



Keep the extra doublers, and divide their effects out by hand. It turns out to 
be simple to reduce from 16 to 4 Dirac fermions, the result being "staggered" 
fermions [T^. I explain this when I discuss Bk below. 



Avoid the Karsten-Smit result by using a random lattice |T^. This is very 



difficult to analyze, because even free fermions must be studied numerically. 
Little progress has been made — see Ref. for a recent study. 

Use "domain-wall fermions" or their descendents. I have no time even to 



introduce the ideas; see Ref. [18| for a review and references to the literature. 
It is controversial at present whether the scheme can be implemented in a 
practical way. 

Most present simulations use Wilson fermions, so I will briefly explain what these 
are. A mechanistic way of understanding doubling is to note that the "forward- 
backward" derivative of naive fermions is small both for a smooth function and for 
one that is smooth except that it alternates in sign. On the other hand, a discrete 
second derivative will be small for the smooth function, yet large for the alternating 
function. Thus Wilson suggested adding to the action a second derivative term 

Sw = J2 7;Qn{Un,^,qn+^, " 2g„ + f/I_^,^g„~^) , (48) 

where r is a parameter. The resulting free propagator is (I leave this as an exercise) 

G{k) = -^^+^-^-f^') , (49) 
s2 + (m+|A;2)2 

where k = 2sin(^). If k^ = tt, then = but k^ = 2, so that the additional 
pole picks up an effective mass rrigff = m + 2r. Thus if one keeps r flnite in the 
continuum limit, the effective mass stays flnite, even if m = rriphysa 0. Thus the 
effective physical mass becomes inflnite. The same applies to all the other doublers. 
Thus we expect them to decouple from the theory in the continuum limit, leaving 
a single Dirac fermion. This slightly sloppy analysis can be conflrmed by studying 
the position of the poles in G{k). 
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In the presence of interactions, the doublers stiU decouple, but they do cause a 
renormahzation of the parameters of the original Lagrangian. This is the standard 
result of the Applequist-Carrazone decoupling theorem. In particular, the quark 
mass m gets additively renormalized. This is not a surprise, since Sw explicitly 
breaks chiral symmetry, so there is nothing special about the value m = 0. Chiral 
symmetry should, however, be restored in the continuum limit, because the Wilson 
term is a discretization of aqd^q, which vanishes when a -— > 0. 

In practice one uses the value r = 1 in simulations. One reason for this is 
that the lattice theory then satisfies reflection positivity, guaranteeing that one can 



construct an hermitian positive Hamiltonian [19 



4 Simulations 

To illustrate what is involved in simulations, consider the pion two-point function 

{uj5d{n)dj5u{p)) = f[dU]Y[[dq][dq]e-^'-^''~^'^'u^5d{n)d-f5u{p) (50) 

- J[dU] n, detip + m,)e-5^Tr Up + m,)-^^,{p + m^);^^. 



/[(it/]n,det(|) + m,)e-^. 

Here Sg is the lattice gauge action, n and p are lattice sites, and p+m is the complete 
lattice Dirac operator appearing in the sum of naive and Wilson terms Sn + Sw 
Simulations are done using the form on the second line, i.e. after integrating out 
the Grassman fields. One is left with the functional integral over gauge fields, with 
measure 

dfi = [dU] n det{p + mq)e~^« . (51) 
g 

To reduce the problem to a finite number of degrees of freedom one uses a lattice 
of finite extent in Euclidean time and space. The finiteness in time actually corre- 
sponds to putting the system in the canonical ensemble at temperature T = l/{Nta). 
For the studies I will consider, Nta is large enough that T ^ Aqcd, and the system 
is effectively at zero temperature. One must choose boundary conditions on the 
fields. To obtain the canonical ensemble these must be periodic (antiperiodic) in 
the time direction for gauge fields (fermions). In the spatial directions any choice 
can be made, though typically periodic boundary conditions are used for all fields. 
One then generates a set of "gauge configurations" (i.e. a U for every link) ac- 
cording to the measure rf/i, and calculates propagators, {p + m)~\ on each of the 
gauge fields. These are joined together into traces such as that in Eq. ^ Finally, 
the resulting correlator is averaged over the "ensemble" of configurations, giving a 
statistical estimate of the desired quantity. 

What is the magnitude of this task? I am involved in calculations on = 32'^ x 64 
lattices. Since each U integral is over the eight dimensional group manifold of 
SU{3), we are attempting to evaluate approximately an 8 x 4(links/site x V ^ 
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7 X 10^ dimensional functional integral. The Dirac operators are complex matrices 
of dimension 

[3(color) X 4(spin) x Vf = [2.5 x 10^]^ . 

They are, however, sparse involving connections only between nearest neighbors. 
Furthermore, if one fixes the site n in Eq. |5^, one need only calculate a single 
column of the inverse (and a single row, but this can be obtained without extra 
work) . In this way one obtains the two point function from n to all other points p. 
When lattice practitioners talk about calculating propagators, they almost always 
are referring to such a truncated calculation. The equation to be solved for the 
propagator G is 

J2{p + m)npGp = Sn , (52) 
p 

where Sn is the "source". This might be 5n,no5 or an extended function. 

Numerical simulations, then, break up into two parts: generating configura- 
tions, and calculating propagators. The latter is done using a variety of standard 
algorithms]^, conjugate gradient and minimal residue typically being preferred for 
staggered and Wilson fermions respectively. Convergence can be improved by pre- 
conditioning. Nevertheless, present algorithms for calculating propagators take a 
number of iterations which is, for small m, roughly proportional to 1/m. This slow- 
ing down is a big problem in practice, as it forces us to work with quark masses 
much heavier than the physical up and down quark masses. Important steps towards 
alleviating this problem have been taken using "multi-grid" algorithms 



There are various methods for generating configurations — Metropolis, Langevin, 
Molecular dynamics, . . . — none of which I have time to discuss in detail. All use 
importance sampling, i.e. they move through the space of possible configurations 
in such a way that the resulting distribution is weighted directly according to the 
measure dfi. If one defines an effective action by 

exp(-^,^(t/)) =ndet(^ + m,)e-^« , (53) 
q 

then, as in any problem with a large number of degrees of freedom, S^jj is very 
highly peaked in U space. It would be utterly hopeless to generate configurations 
uniformly in U space, and then include exp{—Seff) in the integrand. All the methods 
work by moving towards the minimum of Sejj but including some noise to keep the 
appropriate distribution around the minimum. 

Most of the methods will only work if exp(— S'ejff) is positive, in which case it 
can be interpreted as a probability distribution on U space. The gauge action is 
automatically real. In the continuum limit, for a vector representation of fermions 
(which, as we have seen, is what we are forced to use on the lattice, and, in any case, 
is appropriate for QCD), the determinant of the Dirac operator is real and positive 
definite (for m > 0). This is because the eigenvalues come in complex conjugate 
pairs, with eigenvectors related by multiplication by 75. On the lattice, with Wilson 
fermions, the determinant can be shown to be real, but is not necessarily positive. 
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Thus to obtain a positive measure one must simulate degenerate pairs of quarks, in 
which case the measure contains the square of the single quark determinant. 

Almost all algorithms are based on making a change 5U in the gauge fields, and 
then calculating the change in Sejf 

SSeff = - 5] 5Tr \n{p + rUg) + SSg (54) 
q 

= - Jr [6U{p + nig)-']- [3 J2 ^^5 (Re Jrn) . (55) 

q □ 

To calculate the variation of the gauge action involves only a local calculation, i.e. 
one involving links close to that being changed. In contrast, the variation of the 
fermionic part of the effective action involves a propagator and is non-local. Thus 
it is doable (at least one doesn't have to calculate the determinant itself!), but it is 
slow, as it takes a number of operations that grows proportional to V. 

The present algorithm of choice for simulating QCD is the "hybrid Monte Carlo" 



algorithm . One can estimate that the time it takes to generate an independent 
configuration is roughly CPU oc V^^^m~^^^^. As the quark mass decreases, so does 
the pion mass, oc m. Naively, it is necessary for the lattice length to exceed the 
pion Compton wavelength by a factor of a few, so L oc l/rriTr- Using this, one finds 
CPU oc m~^'^^. This is an asymptotic estimate, at best only roughly applicable 
for today's lattice sizes and quark masses. But it shows why it becomes rapidly 
more difficult to simulate lighter quarks. This is why you rarely hear of simulations 
with quarks lighter than m^/S, and why progress is slow, even with an exponential 
growth in computer power. 

There are various ways in which one can overcome this scaling law. First, if 
one uses an improved action, one can get away with a larger lattice spacing (for 
a given size of systematic error), and thus a smaller number of lattice points for 
a fixed physical volume. Second, a lot of the difficulty comes from the slowing 
down of propagator inversions as m — 0, and this should be avoidable. This is the 
aim of the multigrid inversion algorithms. Third, once the quark masses become 
small enough that the pions are light, they interact weakly, and it aught to be 
possible to use chiral perturbation theory to account for errors introduced, say, by 
not increasing the volume as L oc With the presently available 10 Gigafiop 

machines, however, and with the notable exception of studies of QCD at finite 
temperature, systematic calculations of phenomenologically interesting quantities 
are not possible in QCD. 



4.1 Quenched QCD (QQCD) 

To make progress one must make an approximation, and what is used is the so-called 
"quenched" approximation 

djJQCB = [dU] n det(^ + m,)e-^'' — > c?yUQQCD = [dU] exp{-Sg) . (56) 
q 
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In other words, we set the fermion determinant to a constant. This amounts to 
throwing away internal quark loops, while keeping the valence quarks, which now 
propagate through a modified distribution of gauge configurations. For this reason, 
it is sometimes called the "valence" approximation. Using QQCD reduces CPU 
requirements by a factor of ~ 10^ — 10^ with present values of a and m. Furthermore, 
the time to generate new configurations only grows as CPU oc V^/^ oc m so that 
it is easier to go to smaller quark masses. 

Throwing away quark loops is a drastic approximation, the effect of which I 
discuss in more detail below. Nevertheless, it is an approximation certainly worth 
studying, because it retains the essential non-perturbative features of QCD, confine- 
ment and chiral symmetry breaking. A quark propagating through the ensemble of 
quenched gauge fields picks up an effective mass, just as it does in QCD, and binds 
to form hadrons. The effective mass and the details of the binding will differ from 
those in QCD, but, given the success of the quark model in describing the hadron 
spectrum, one would expect the quenched spectrum to be qualitatively similar. 

How well do we expect the QQCD to reproduce the spectrum of QCD? Since we 
are stuck with QQCD for a while to come, it is important to try and estimate such 
quenching errors. Let me begin with a rough estimate. One of the unphysical effects 
of quenching is that resonances in QCD, e.g. the p meson, become stable states in 
QQCD. This is because internal quark loops are necessary to obtain the on-shell 
intermediate states (e.g. tttt in the case of the p) which give rise to the imaginary 
parts of the propagators, and thus to the width of resonances. Discarding these 
intermediate states, however, affects not only the imaginary part, but also the real 
part of the propagator. In other words, not only is the width of the state changed 
(to zero) but also the mass is shifted. The most naive estimate is that 6m ~ 5r = P. 
This mass shift will not be uniform in sign or magnitude, since it depends on the 
available thresholds, possible cancellations, etc. It may, in fact, be small for the p 



22l. But it wiU, I expect, distort the spectrum at the 10% level (lOOMeV/lGeV) 
in general. 

Can one make a better estimate of the effects of quenching? In particular cases, 
I think one can. One example is quarkonia (cc, 66), which I touch on below. Here 
I describe a somewhat more systematic method applicable to the properties of the 
pseudo-Goldstone bosons (PGBs), the vr's, K^s and rj. 



The starting point is the formulation of QQCD introduced by Morel[23 



Zqqcd = J [dU] [dq] [dq] [dq] [t;|]e-^«e^(^+™)«e-^(^+"^)^' . (57) 

Here g is a ghost field: a commuting spin-| variable. I have written the partition 
function for only one fiavor — in general there is a ghost degenerate with each quark. 
Equation |5^ works because the ghost integration yields an inverse determinant 
which cancels that from the quark integration. In other words, internal ghost loops 
cancel internal quark loops exactly. This formulation shows why QQCD is a sick 
theory — ghosts have the wrong connection between spin and statistics, and thus, 
in Minkowski space, there will be violations of causality. One cannot avoid the 
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Figure 6: Contributions to the r( propagator. Lines are quark propagators. 

problems by considering correlation functions of external states composed of quarks 
alone, because particles containing ghosts will appear in intermediate states. Note, 
however, that, even though the Minkowski theory may be sick, the Euclidean version 
is well defined. This only requires that the eigenvalues of ^ + m have positive real 
parts, so that the bosonic functional integral converges. This is an example of a 
Euclidean theory which is not reflection positive. 

Morel's formulation is the starting point for the development of "quenched chiral 
perturbation theory" (QChPT) [||, Because of the ghosts fields, QQCD has a 
larger chiral symmetry than QCD, namely 



Here S't/(3|3) is the graded group of special unitary transformations of three com- 
muting and three anticommuting objects. Numerical evidence suggests that this 
group is spontaneously broken down to its vector subgroup, yielding not only the 
usual pseudo-Goldstone bosons (PGBs), but also partners in which the one or both 
of the quark and antiquark fields are replaced by ghosts. One can develop an effec- 
tive Lagrangian for this theory, analogous to the usual chiral Lagrangian. Using this 
one can calculate the effect of loops of PGBs, which give rise to non-analytic terms 
generically referred to as "chiral logs". These logs are, in most cases, different in 
QCD and QQCD, and one can use the difference to estimate the effect of quenching. 
In particular, it should be possible to give a rough ordering of quantities according 
to the size of quenching errors. I give some examples below. 

There is an important qualitative difference between QChPT and ChPT (chiral 
perturbation theory for QCD). In QCD, the Vj is not a PGB, since the would-be 
f/(l)^ symmetry is anomalous. The Vj gets additional mass from diagrams in which 
the quark-antiquark pair annihilate into gluonic intermediate states. This "hairpin 
vertex" gets iterated to give rise to the mass term, as shown schematically in Fig. 
|[ But in QQCD, all terms except the first two are absent, and so the r]' remains 
light, effectively a PGB. The would-be mass term becomes a two-point vertex. The 
77' must be included in the quenched chiral Lagrangian, along with this additional 
vertex. The quenched chiral Lagrangian of Ref. provides an explicit realization 
of these rather handwaving diagrammatic arguments. As explained below, the 
existence of a light Vj leads to a number of unphysical effects. In particular, there 
is an Vj cloud surrounding all hadrons containing at least one light quark. 



(58) 
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5 Numerical Results from quenched QCD 



In this section I present the evidence which shows that QQCD does give a reasonable 
approximation to the real world. This is essential if we are to have any hope of 
using QQCD to calculate phenomenologically interesting matrix elements. I will 
also discuss in more detail the size of quenching errors. 

5.1 Confinement 

There is a simple criterion for confinement in QQCD: evaluate the potential energy 
V{R) of an infinitely heavy quark [Q) and antiquark [Q) as a function of the distance 
R between them. The standard picture of confinement has a tube of color electric 
flux joining the quark and antiquark, a tube which simply elongates as the pair are 
separated. This suggests that V{R) oc R for distances significantly larger than the 
width of the tube. A better way of stating this is that the force, F = —dV/dR, is 
expected to asymptote to a constant at large R. The magnitude of this constant is 
called the string tension, k. Using the force avoids the problem that V{R) contains 
the self energies of the quark and antiquark, which are R independent, but divergent 
in the continuum limit. 

This criterion for confinement is not useful in QCD, because the flux tube can 
be broken by the creation of a qq pair, leaving a Qq and qQ meson having an energy 
independent of R. This process does not occur in QQCD because it requires internal 
quark loops. 

To evaluate V{R) one proceeds as follows. Create the QQ pair at time r = 
using the gauge invariant operator 

Q{R)L{R,0)Q{0) , (59) 

The ordered integral in L (Eq. ^) can follow any path, or one can average over a 
number of paths so as to maximize the overlap of the operator with the QQ state. 
Next, destroy the pair at a later time r = T using the conjugate operator 

Q(0)L(0, R)Q{R) = Q{0)LiR, 0)tQ(i?) . (60) 

In this way one has constructed a correlator which, for large T, should fall as 
exp(— V(i?)T), where V{R) is the energy of the lightest state of the Q + Q + glue 
system. 

To evaluate the correlator we need the heavy quark propagator. In Minkowski 
space, an infinitely heavy quark just maintains its velocity, so if it starts at rest, as 
here, then it remains so. That is why the final operator in Eq. ^ is chosen with 
the Q field at the same site as the original Q. The current of a static quark does 
have a time component, which couples to Aq, resulting in a phase for the propagator 
Pexp[—ig J dtAolt)). There is also the kinematical phase exp(— zMgt), but this does 
not depend on R and thus does not contribute to the force, and can be dropped. Ro- 
tating to Euclidean space, the line integral remains a phase: P exp{—ig J drA^lr)). 
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Figure 7: The heavy quark potential, in lattice units. The short distance points have 
been corrected for lattice artifacts using the lattice Coulomb propagator. 



For the antiquark, the line integral is in the opposite direction. Combining these 
propagators with the line integrals from Eqs. ^ and we obtain a Wilson loop, 
W. It is roughly rectangular, having straight segments in the time direction, while 
the spatial paths (determined by the choice of L{0,R)) can wiggle, though they 
must remain in a single time-slice. Putting this all together, we expect 

(^W) Ce-^(^)^ ^ Ce-'^^^ , (61) 

where C is a constant. The last result is the well-known "area-law" : the expectation 
value of a rectangular R x T loop falls exponentially with the area of the loop if 
there is confinement. 

It was established long ago, using the strong coupling expansion, that the area 
law holds for large g"^. This leaves the question of whether the result extends to 
the continuum limit, — * 0, i.e. whether there is a phase transition at finite g^ 
which breaks the analytic connection between weak and strong coupling. Numerical 
evidence to date suggests that this does not happen. For example, I show in Fig. 
the results for the potential obtained by Bali and Schilling p^] at /5 = 6.4, which 



corresponds to a ^ 0.06 fm. The horizontal scale is in units of a, and thus extends 
well beyond Ifm. The linear behavior of V is clear, starting from ~ 0.5 fm. A 
pleasing feature of this result is that one can see, in the same calculation, both the 
long distance, non-perturbative physics of confinement, and the short distance per- 
turbative Coulomb potential. This is actually necessary if one is to calculate weak 
matrix elements, for one must match onto the continuum using perturbation theory 
at short distances, while simultaneously including the long distance contributions. 
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5.2 Charmonium and Bottomonium Spectra, and extracting 

as 

I now turn to the spectrum. After integrating out the fermions, and making the 
quenched approximation, a general meson two-point correlator becomes 




time 



X 



where the "blob" at r = represents an initial extended source, made gauge invari- 
ant in some way, and the blob at r represents a similar "sink" . The lines joining the 
blobs are quark propagators in the background quenched gauge field. If the meson 
is a fiavor singlet (e.g. cc), then there is a second diagram in which the quark and 
antiquark annihilate into intermediate gluons. For heavy quark systems these are 
suppressed by powers of as{mq), and can be ignored. The same is not true for light 
quark systems; in particular, the r]' mass comes from such diagrams. 

I first discuss the cc system. This has been studied on the lattice with great care 
by the Fermilab group [^. Compared to the light quark spectrum, the cc system 
has several advantages 

• cc states are smaller than light quark hadrons, so one can use lattices with 
smaller size in physical units. 

• The CPU time needed to calculate c-quark propagators is less than for light 
quarks. Since, in QQCD, calculating propagators consumes most of the com- 
puter time, this allows one to use more lattices, and thus reduce statistical 
errors. These errors are further decreased by that fact that the intrinsic fiuctu- 
ations of c-quark correlators from configuration to configuration are typically 
smaller than for light quarks. 

• The cc system is reasonably well described by potential models, so one has a 
way in which to estimate systematic errors, in particular that due to quenching. 

A potential disadvantage is that the charm mass in lattice units is not small, e.g. 
rrica ~ 0.75 at a = 0.1 fm [jS ~ 6). This might lead to significant discretization 
errors, proportional to powers of m^a, and thus require the use of a smaller lattice 



spacing. This turns out not to be true, as long as one uses an improved action p8 |. 
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Charmonium is thus a system where aU of the lattice errors can be studied, and 
to a large extent controlled. In particular, finite volume and lattice spacing errors 
are small. It turns out, as I discuss below, that similar control is possible for the bb 
system, using non-relativistic QCD (NRQCD) to describe the b quark. The onia are 
thus a good choice for accurately measuring the lattice spacing. As I will explain, 
this can be turned into a prediction for a^. 

The lattice spacing is obtained by comparing a quantity measured in lattice units 
to its physical value. An excellent choice for the physical quantity is the splitting 
between the spin-averaged IP and IS" levels in onio. This is because the splitting 
is almost the same in cc and bb systems (457 and 452 MeV, respectively), so that 
one need not worry if the lattice quark masses differ slightly from the experimental 
values. To extract a one uses 

where (awiis',p)iat are the dimensionless masses one obtains from the lattice simula- 
tion. The absence of corrections of O(a^) in this relation is a convention. If we used 
other physical quantities we would obtain values of a differing by such corrections. 

This method gives accurate values of a for various values of (7^, or equivalently 
as a function of a. From the point of view of perturbation theory, the lattice is just 
an ultra-violet regulator, albeit a messy, rotationally non-invariant one. Roughly 
speaking, it restricts momenta to satisfy \p\ < vr/a. Thus it can be related to 
couplings defined in other schemes, for small enough a. For example, it is related 
to the MS coupling 

g\a) = gl^{fi = 'K/a) [l - 0.31g' + 0{g')] . (63) 

(Neither the scheme nor the scale of the g^ in the correction is determined at this 
order.) The idea is to use this result to obtain gj^{7i/a), and then use the beta- 
function to run this coupling to a standard scale, for example mz- 

An important test of the calculation is that one should obtain the same result 
for gy^{mz) starting at different values of g^. The extent to which this is not true 
is an indication of the size of errors due to truncating perturbation theory in Eq. 
|63i from truncating the beta- function when running to m^, and from discretization 
errors. It turns that the combined effect of these errors is smaller than the statistical 



errors 
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The calculation cannot be carried out exactly as just explained. The g term in 
Eq. ^ is ~ 30%, since ~ 1 in lattice calculations. Thus there are likely to be 
large corrections to Eq. coming from unknown higher order terms. How, then. 



do we convert the well calibrated lattice into a result for a^jg? We need a non- 
perturbative way of obtaining «r?o. Such a method has been suggested by Lepage 



and Mackenzie [11 



There are several parts to their method. First, we recall from our experience 
with perturbative QCD that a^Asif^) ^ reasonable expansion parameter, as long 
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as one chooses a scale /i appropriate to the process under consideration. Equation 
|63| then imphes that aiatia) = 5'^(ct)/4vr will be a poor expansion parameter. This 
expectation is borne out by numerical results. For example, ratios of small Wilson 
loops are not well represented by either first or second order lattice perturbation 
theory — the 0{g^) term is roughly half the size of the needed correction at = 1. 
If, however, one expands these quantities in terms of a^s, and chooses the scale 



according to a prescription explained in Ref. [|Tl[ , then the leading order term does 
much better (because ajj^ > aiat), and the second order expression works very 
well. See Ref. [|lll for other examples. The moral is that perturbation theory 
for short distance lattice quantities is in good shape, as long as one chooses the 
correct expansion parameter. It is also noteworthy that Lepage and Mackenzie 
have understood the source of the large correction in Eq. it comes from extra 
"tadpole" diagrams that occur in lattice perturbation theory. 

I have skipped over an important detail in the previous paragraph. Lattice 
perturbation theory is rapidly convergent for almost all quantities when expressed 
in terms of a^^ at an appropriate scale. But how to we determine a^jg, given the 
need for higher order terms in the relation Eq. |6^? Lepage and Mackenzie suggest a 
non-perturbative definition in terms of the numerical result for the average plaquette 

. _iMI[H).|:,i,o(/)) (64) 

1 1 

0.37 + 0(a). (65) 



a|;jg(3.4l7r/a) ap 

In other words, define an auxiliary coupling constant by the first equation, and 
relate it to aj^ using the perturbative result in the second line. Note that the 
correction term in the second relation is small, since 1/a ~ 5 — 10. The scale 
SAln/a takes into account a subset of the two-loop corrections. It is using ajj^ 
defined in this way that Lepage and Mackenzie find lattice perturbation theory to 
be well behaved. In effect this method uses the numerical data itself to sum the 
leading tadpole diagrams to all orders. 

In summary, using the result for the lattice spacing from Eq. B2|, together with 



the result for the average plaquette, Eqs. |6J and ^ give us at the known scale 
q = SAln/a. We can then run the result up to niz- 

By far the largest error in the result is that due to the use of the quenched 
approximation. The Fermilab group has developed a method to estimate this error, 
based on the fact that the cc system is well described by a potential model. The 
potential in QCD differs from that in QQCD, but we have a reasonable idea of the 
form of this difference, so we can subtract its effects. It is useful to define a coupling 
constant in terms of the potential in momentum space 

V{q) ^ -Cpa^v'\q)/q' ■ (66) 

The superscript gives the number of fiavors in internal quark loops. At short dis- 
tances this behaves like any other coupling constant, and in fact is close to that in 
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the MS scheme 

— rrr^ = rrr^ - 0-822 + 0{a) . (67) 

At long distances it must lead to a confining potential. A useful form for a which 
interpolates between these two limits has been given by Richardson. Now, since 
one adjusts the scale so that the IP — IS splitting matches experiment, it must 
be that QCD (uf = 3) and QQCD {uf = 0) potentials are similar at the scale of 

typical momenta in low lying cc states, q* = 0.35 — O.TGeV, i.e. av\g*) ~ <^v\l*)- 
But the two couplings run differently with q. In particular, for large enough g, 
the QQCD coupling decreases more rapidly because there is no fermionic screening. 
This means that ay^(7r/a) — ay\'K/a) > (assuming rr/a > q*). The difference is 
calculable given an assumed form for the potential. Using Eq. |6^ we can convert 
this to a difference between a^{TT/a) and a^(7r/a). The former we have already 
determined, so we obtain the latter, which we then run up to mz- 

The resulting correction is substantial — a 26% increase in a(^(5GeV) for char- 
monium. This correction itself is uncertain, due to uncertainties in the matching 
scale g* and in the chosen form of the potential. The resulting uncertainty in 05 is 
estimated to be ~ 8%, slightly larger than that coming from the neglect of higher 
order terms in the perturbative relation between ap and a^jg. A clear discussion of 



all these issues is given in Ref. p8 



The net result is 28 



ag^(m^) = 0.110 ±0.008 . (6^ 



A similar analysis using NRQCD to study the bb system yields a consistent result] 
0.112 ± 0.004. The fact that the final results agree is a test of the method of 
estimating quenching errors, which are different in the two onia. It is appropriate 
that the result be included in the latest Review of Particle Properties fS^^ . 



Recently, the first results including quark loops have been obtained, using two 
moderately light fiavors of quarks in the loops P^, PT| . For a review see Rev. ||33|| . 
Results for three light fiavors can now be obtained both using the methods just 
described, or simply by extrapolating from Uf = and Uf = 2. The two methods 
yield consistent results, with the latter giving the smaller errors: [ 



a^^imz) = 0.115 ± 0.002 . (69) 

This is a very impressive result, which is consistent with the latest world average 
experimental value a^(m^) = 0.117 ± 0.005. My only concern is whether the 



error fully accounts for the uncertainties in the extrapolation of ruu and to their 
physical values. 

5.3 Light hadron spectrum 

I want to begin my discussion of the spectrum of hadrons composed of u, d and s 
quarks by clarifying how, were we able to simulate QCD, we would determined the 
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correct lattice quark masses and extract a. This must be done carefully because 
the measure of the functional integral depends, through the Dirac determinant, on 
the quark masses. For simplicity of presentation, I assume that m„ = md = mi. I 
will also assume that we have picked a value of (3 = G/g"^ large enough that 0(a) 
errors can be ignored. One would remove such errors in practice by simulating at a 
number of values of f3 and extrapolating to a = 0. 

Having chosen (3, we then calculate numerically three lattice masses, e.g. 

amprotonl'^'^/; ct^^s) ; amT^{ami , aiTLs) , and am xiaTrii, arris) (70) 

as a function of the lattice masses m;a and m^a. We adjust these lattice masses 
until the two mass ratios agree with experiment 

amT,{W7}4,(ifn^) 135 amK{oMi,aWr^) 498 ^^^^ 
am^^oton{aMl,am;) 938 ' amp,oton{am,afrT;) 938 ' 

Hopefully these equations have a solution! With the lattice masses fixed we can 
now extract the lattice spacing, using 

amproton{oMi,Wm^) = 938MeV . (72) 

Any other quantity (mass, decay constant, . . . ) that we calculate can now be 
predicted in physical units. Clearly, any three dimensionful quantities can be used 
to carry out this program. 

Knowing a also allows us to extract the quark masses in physical units. These 
are quark masses renormalized in the lattice scheme at the scale a. They are per- 
turbatively related to more familiar masses, such as those in the MS scheme. 

The same procedure applies to QQCD, but is much simpler to implement. This 
is because the measure is independent of the quark masses, so one need only generate 
a single ensemble for each (3. The flip-side of this simplicity is, of course, that one 
is throwing away important aspects of the physics. 

Unfortunately, it will be some time before this procedure can be followed for 
QCD. For the moment a systematic study of the spectrum is restricted to QQCD. 
The most thorough analysis is that of the IBM group, using the GF-11 computer 
built by Weingarten and collaborators 0. It sustains ~ 7 GFlops for these calcula- 
tions. I will present an outline of their results, with my major focus being the issue 
of the reliability of the quenched approximation. 

The main features of the study are 

• Masses are calculated for light hadrons (those with the quantum numbers of 
the TT, p, N and A), using degenerate quarks. A range of quark masses is 
used, extending down to about mga/S, where is the physical strange quark 
mass. The results are then extrapolated to the physical up and down quark 
masses. An example is shown in Fig. ^. I return to the reliability of these 
extrapolations below. 
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The use of degenerate quarks means that a direct calculation of the strange 
meson and baryon masses (except for the VL~) is not possible. There is no 
fundamental obstacle to doing such a calculation; after all, quark propagators 
with different masses have been calculated. The practical problem is just that 
the number of non-degenerate combinations becomes large. 

The IBM group make predictions for the strange hadrons based upon the 
assumption (well supported by the experimental data itself) that, in QCD, 
the masses of mesons and baryons (squared masses for PGBs) depend linearly 
on the masses of the valence quarks of which they are composed. This leads, 
for example, to the result 

ms + — TT^N ~ mN^rris) , (73) 

where the quantity on the right hand side is the mass the nucleon would have if 
= = rris for the valence quarks. This latter quantity is easily calculated 
in QQCD. 

Finite volume effects are studied by comparing results on lattices of different 
sizes. Previous studies have shown that for La < 1 — 1.5 fm there are significant 
finite volume corrections . All the lattices in the IBM study are larger than 



this, so the finite volume corrections turn out to be small. 

The extrapolation to a = is done using three lattices, having roughly the 
same physical volume (L ~ 2.4fm), but different lattice spacings {(3 = 5.7, 
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Figure 9: Extrapolation to a = (for sinks 0,1, and 2 combined). All masses have 
been previously extrapolated to the continuum limit. Dots at a = are observed 
values. 

5.93 and 6.17, corresponding to a ~ 0.15, 0.1 and 0.07fm). Figure |^ shows 
an example for rri]^ /nip and "(ms + m= — mj^)/mp" (obtained using Eq. [731) . 
Since Wilson fermions are used, discretization errors are 0(a), and linear 
extrapolation is assumed. This plot does not include a small shift due to 
finite volume corrections. 

• Statistical errors are reduced using large ensembles of lattices, typically 200. In 
addition the signal is improved by creating the hadrons using extended sources. 
Figure ^ shows about the worst signal. What is important is that the effective 
mass reaches a plateau before it disappears into noise. The fitting range is 
chosen by an automatic procedure, and is shown by the dotted vertical lines. 
The data sample is large enough that fits using the full correlation matrix are 
stable, and errors are estimated using the bootstrap procedure. The resulting 
mass value is shown by the solid horizontal line, and has statistical errors of 
about 2%. 

The final results are shown in Fig. 10. The predictions all agree with the experimen- 
tal numbers to within 6%, or less than 1.6 standard deviations. QQCD seems to 
work extremely well! Indeed Ref. ^ interprets this success as indicating both that 
QCD is the theory of the strong interactions, and that the quenched approximation 
works well, at least for masses. 
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Figure 10: Final results for the light hadron spectrum. Masses are measured in units 
of rrip. Experimental numbers are denoted by crosses. 



5.4 Quenching errors in the Ught hadron spectrum 

The IBM study is an impressive piece of work, one which sets the standard for 
future simulations. Only a few other quantities (a^, nib and Bk) have been studied 
so thoroughly. I am not, however, fully convinced of the conclusions, i.e. that the 
spectra of QQCD and QCD differ by less than 6%. As I will explain, there are rea- 
sons to expect typical deviations to be 10-20%. There are then three possibilities: 
(i) the reasons I will present are not valid; (ii) the reasons are valid, but the devi- 
ations turn out fortuitously to be smaller for the light hadron spectrum; and (iii) 
the assumptions made in order to do the various extrapolations in the IBM study 
are not valid, and the actual answers differ more substantially from the physical 
spectrum. I am biased, so I expect the answer to be a combination of (ii) and (iii), 
but only further simulations will resolve the issue. 

The two extrapolations which might need further refinement are those in lattice 
spacing and in quark mass. As for the former, one expects corrections of the form 1 + 
aAi-|-(aA2)^ + . . ., where Ai and A2 are non-perturbative scales. Linear extrapolation 
has been assumed (A2 = 0). But if A2 ~ nip, then the quadratic term would be 
a 25% correction at the largest values of a used in the extrapolations, and could 
not be ignored, illustrated by For example, for the nucleon data in Fig. ^ such a 
quadratic term could lead to an extrapolated value m^v/mp ^ 1.35, in contrast to 
the result, 1.28, from a linear fit. Such large values for A are not unreasonable — they 



have been seen in the calculation of with staggered fermions . 

The chiral extrapolations, such as those of Fig. H, have been done linearly using 
the lightest three mass points. Looking at the nucleon data, there is some evidence 
of negative curvature, as has been also seen in other simulations. As explained 
below, we expect there to be an term with a negative coefficient. Including such 
a term would lead to a slightly lower extrapolated value. 
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Figure 11: Quark diagrams leading to pion and clouds around mesons. 

My point here is not that the extrapolations are necessarily wrong, but that 
there are reasons to examine them more carefuUy, in particular by accumulating 
data at more values of a, and at smaller quark masses. 

Let me now explain the reasons why I expect the spectra of QQCD and QCD to 
differ by 10-20%. The essential point is that hadrons in QQCD have very different 
"pion" (vr, t] and f]') clouds than in QCD. To understand this, consider Fig. 
pT| which shows "quark-flow" diagrams for processes leading to pion loops. The 
justification for the use of such diagrams is discussed in Refs. [2^, 2^. The first 



diagram is present in QCD, but absent in QQCD. The second is present in both 
theories, but in QCD it involves the fiavor singlet r)', which is heavy and thus does 
not contribute significantly to the long-distance cloud. In QQCD, by contrast, the 
T]' remains light, as discussed above, and this diagram leads to the appearance of 
an T]' cloud around hadrons. So what happens is that the pion cloud surrounding 
mesons in QCD is replaced by an rj' cloud around quenched. These two clouds are 
not related — the rj' cloud is simply a quenched artifact. To the extent that particle 
properties in QCD depend upon the composition of the cloud, they will be altered 
in QQCD. 

The extent of the alteration can be investigated using (quenched) chiral per- 
turbation theory. One of the most striking results is that the chiral limit is more 
singular in the quenched theory. For example, the QChPT result for the pion mass 
is 



24, 25 



,2 



2fxmg l-2(51n(m^/A)+0(m^ln(m^/A')) . (74) 



Here 5 is a constant proportional to the rj' two point vertex in the right-hand diagram 



of Fig. |TT]. The term proportional to 5 diverges as m^r 0. This is contrast to the 
corrections in QCD, which are proportional to m^ln(m7r), and vanish in the chiral 
limit. Thus the r]' loops of QQCD introduce a sickness that is entirely an artifact. 
Indeed, the chiral expansion ceases to make sense once gets small enough that 
the "5 term" is comparable to l.| In fact, a general feature of QChPT is that there 
are corrections which diverge as one or more quark masses are sent to zero. This 
suggests that one cannot hope to use QQCD below a certain quark mass. This 



§For a pion made of degenerate quarks, one can resum the leading terms, [g5[ |3^ i.e. those 
proportional to (Shim^)^ . It is not known how to do this in general. 
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Figure 12: Searching for the t]' cloud. All masses are in lattice units. The two fits 
lead to the values of 6 shown. 



"critical" mass will likely depend on the quantity being calculated. 

To estimate the critical mass one need to know 6. In QCD, the two-point vertex 
leads to the major part of m^/, the remainder due to quark masses. Assuming the 
vertex is the same in QQCD as in QCD, one finds 6 ~ 0.2. In principle, one need 
not appeal to QCD, but rather can calculate the rj' two-point function in QQCD 
and directly extract 6. This is a difficult quantity to measure, requiring propagators 



from many sources. A calculation has been recently completed fS^], however, albeit 
at a rather large lattice spacing (a ~ 0.15 fm). The result, 6 ~ 0.15, is close to the 
QCD-based estimate. The smallness of 6 justifies the use of perturbation theory, 
and implies that critical masses turn out to be small. 

There is now some numerical support for QChPT. Figure |12| shows a test of Eq. 



f4[ I have taken this from the recent review of Gupta [^^, which provides more 
details of the fits. QChPT predicts that the ratio m^/m^ should diverge at small 
quark masses. This can be tested using staggered fermions, for which the quark 
mass is mult ipli cat ively renormalized, and one knows where nig = 0. The data at 
small quark masses are from Ref. , and come from three lattice volumes. There 
is a pronounced finite volume effect betwee 16'^ and 24^ lattices, but no such effect 
from 24^ to 32'^. Thus the 32^ data is close to the infinite volume limit. A technical 
point, which is crucial numerically, is that the mass appearing in the logarithm in 
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Eq. is that of a pion which is not an exact PGB on the lattice. The best fit 
(soUd hne) gives a value 6 ~ 0.16 — other reasonable fits give different values, but 
all are definitely non-zero. 

A less striking, but equally important, test comes from decay constants (/tt, fx 
etc.) Thes are also expected to diverge for mesons composed of non-degenerate 
mesons when one of the quark masses vanishes. The world's data for such decay 
constants is consistent with the expected QChPT form, if 5 = 0.10(3) ||3^, consistent 
with the two other determinations. To be complete, I should note that there is as 
yet no evidence for a divergence in rn^/iriq for Wilson fermions [0. It is, however. 



much more difficult to do the fit, since the quark mass is additively renormalized, 
and one does not know a priori where = 0. 

Thus I have some confidence that QChPT makes sense, and is applicable to 
present simulations. What does it imply for baryon masses? In QCD the chiral 
expansion of baryon masses is ttin = mo + cml + c'm^ + 0{'m'^) , where represents 
any of the PGB masses, and c, c' are constants. The non-analytic oc m^/^ term 
is due to loops of PGBs, so-called "chiral loops" . Its coefficient is known in terms 
of the PGB couplings to the nucleon, e.g. q-knn- The analytic corrections of O(m^), 
by contrast, involve additional, unknown, parameters. But in the chiral limit, the 
non-analytic term is enhanced relative to the analytic term by one power of m^, so 
one has some predictive power. This is advantageous compared to mesons masses, 
for which the chiral logs of size In m^r are enhanced only by a logarithm over the 
analytic terms of O(m^). 

In QQCD it turns out that for baryons, unlike for mesons, some, though not 
all, of the chiral loops involving PGBs remain pT|. Thus there is still an term. 



though with a different coefficient. There are also new terms, due to r/' loops, 
which are proportional to 5 x m^r, and thus enhanced in the chiral limit. These 
are the analogs of the divergent terms in m^, and are pure artifacts. Jim Labrenz 
and I have calculated the chiral expansion for octet and decuplet baryon masses in 
QChPT|^. To give an example of the results, I use the QCD values for Qt^nn and 



other constants, and include only intermediate octets. We then find (all masses in 
GeV) 

m^r = mo — 0.35(5/0. 15)m^ + 3. 4m^ — 1.5m^ + 0(m^ In m^r) . (75) 

The numerical values are only to be used as rough guides, since the coupling con- 
stants in QQCD will not be the same as in QCD. The most important points are 
qualitative: when plotting vs. there should be curvature at larger masses 
due to the term, and there should be a peculiar behavior at small masses due 
to the term, 



These results are illustrated in Fig. |T3|. This shows the nucleon mass from Ref. 
P at /? = 5.93. I have done three types of fits: a fit of rriQ + bml to the lightest three 
points; and fits of Eq. |7^ to all four points with 6 fixed to be and 0.15, but with 
all three other coefficients free. The results of the fits are given on the plot. All fits 
are reasonable, but the best are the two with curvature. This is qualitative support 
of an term, but is by no means definitive. The observation of curvature depends 



37 




on heaviest mass point, and this is at sufficiently large that higher order terms 
in the chiral expansion could also be significant and give curvature [^^. The data 
provides no evidence for or against an m,r term — as the Figure shows, the "hook" 
which appears for 6 = 0.15 is too small to be important unless one goes to very 
small quark masses. 

In this instance, the variation in mo (the intercept) is quite small: the linear fit 
to ml gives mo = 0.92 GeV, while the fits with 6 = (0,0.15) give (0.89, 0.95) GeV. 
Roughly, then, adding the ml term reduces mAr by 3%, while the adding the m^^ 
term increases mN hy 7%. These are small effects, but they are important given that 
they are larger than the typical statistical errors. It is not clear what to do about 
the m^ term. Since it is an artifact of quenching, it may be best to extrapolate 
ignoring it. But one should certainly include an m^ term. 

I now return to the reason for this digression into QChPT. We want to make an 
estimate of the effects of quenching on baryon masses. We do this assuming that 
QCD and QQCD differ only because the contributions of chiral-loops are different. 
In particular, we assume that all coupling constants in QChPT are the same as in 
ChPT. We apply this method to ratios of baryon masses, and to m^/ f-K- Considering 
ratios removes changes in overall scale between QQCD and QCD. We find that there 
are 10-30% differences between ratios in the two theories ||41||. This is certainly hand- 



waving — the coupling constants could conspire to make the two theories agree more 



38 



closely on the various ratios. But it indicates the magnitude of the typical effect 
due to the difference in the physical composition of hadrons in QQCD and QCD. 
It is because of this general argument that I expect the final quenched spectrum to 
differ more substantially than 6%. We are in the process of extending this analysis 
to other data sets and to the decuplet baryons. 

m. 



For more extensive reviews of QChPT, see Refs. [p8 



6 Anatomy of a calculation: Bk 

1 now turn to applications of QQCD where we do not know the experimental results 
in advance. We have already seen one example — as; from now on I will focus on 
matrix elements of the electroweak effective Hamiltonian. These so-called "weak 
matrix elements" govern weak decays and transition amplitudes. I begin with a 
detailed discussion of the calculation of Bk, not only because it is dear to my heart, 
but also because it clearly illustrates all aspects of such calculations. 

Bk arises when calculating CP- violation in K — K mixing, which is parameter- 
ized experimentally by e. This mixing is caused by box diagrams such as that in 
Fig. |1]. Using the renormalization group (RG), we integrate out the top quark, Z 
and W bosons, and then bottom quark, as well as gluons with momenta exceeding 
the renormalization scale fi. We lower fi down to a scale at which we can match 
onto the lattice calculation, /i ^ n/a ~ 5 — 10 GeV. At this stage the AS = 2, 
CP-violating part of the effective Hamiltonian is, to good approximation 

HM (X Glim [V,IV,1] cifi) [s^,{l+^,)d s-f,{l+^,)d] . (76) 

Here c(/i) is a perturbatively calculable coefficient function, known at present to 
two-loop order. The non-perturbative part of the problem is the evaluation of 
{K\Ti.cs{fi')\K), which is parameterized by Bk{^j) (Eq. |I|). To evaluate this, we 
switch to the lattice renormalization scheme, matching the continuum four-fermion 
operator with a corresponding lattice operator. We then evaluate the matrix el- 
ement of the lattice operator using the numerical methods of lattice QCD, and 
finally convert the result back into one for Bk{i^)- Combined with c(/i) and the 
other constants one obtains an expression for e, from which one can extract the 
value of lm[V,%l]. 

I would like to stress a point which occasionally gets overlooked. What the 
lattice calculation gives us, once we match back onto the continuum, is -Bi^(yu) 
in a continuum scheme of our choice (e.g. naive dimensional regularization) at a 
scale /i which should be near to n/a. We can choose /i to be a standard scale, 
e.g. 2 GeV. This result contains all sorts of lattice related errors, discussed in 
detail below. For phenomenological applications, we must combine it with c(/i), the 
result of a perturbative calculation. This quantity has errors due to truncation of 
the RG equations, and due to the uncertainty in the value of as (or equivalently 
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the value of Aqcd)- These errors have nothing to do with the lattice calculation.|^ 
Since c(/i) oc a~^^'^^{jj,) at leading order, it has become standard to quote results for 
13 X = a7^'^^^(/i)-Bx(A^)- I do not like this practice, for various reasons. The most 
important is that it mixes up errors from different sources. In addition, Bk is only 
fi independent at leading order, so it is not what one actually uses in a two-loop 
phenomenological analysis. And, finally, using 13k amounts to running Bk down to 
a very low scale, that at which a = 1, where I, at least, have little intuition for the 
physics. I propose, instead, that we quote Bk{ij) for a standard scheme and scale, 
just as we do for 0:5. Other methods of calculation (large N^., QCD sum rules, . . . ) 
give Bk at different scales, and will have to be run to the standard scale. If the 
change in scale is small, however, the uncertainties introduced by the running will 
also be small. In this way we can make comparisons between models without the 
overall common errors in c(/i). To do phenomenology, we can take c(/i) from the 
one of the standard RG analyses. 

With that off my chest, let me return to the issue of how we calculate Bk, and in 
particular, how we estimate and reduce the errors. The sources of errors are these. 

• Numerical method. Statistical errors are now at the 1-2% level. 

• Matching continuum and lattice operators. With staggered fermions, the er- 
rors from neglecting two- and higher loop terms in the matching are small, 
1-2%. 

• Making sure the result has the correct behavior in the limit mK 0. This is 
much simpler using staggered than Wilson fermions. 

• Extrapolating to the physical kaon (containing a highly non-degenerate s and 
d) from the lattice kaon (with degenerate, or nearly degenerate, quarks of 
mass ~ ms/2). 

• Finite volume effects. It turns out that these can both be estimated theoret- 
ically to be very small (< 0.5%), and are numerically observed to be smaller 
than the statistical errors on lattices of size 1.6 — 2.4fm across. 

• Errors due to the use of the quenched approximation. 

• Extrapolating to a = — possibly using improved actions. 

I discuss the most important of these in turn. 

^There is a small correlation in the errors, since the lattice-to-continuuni matching coefficients 
depend on as, but this is a minor effect. 
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6.1 Numerical method 



Staggered fermions give the most accurate numerical results for Bk, and I will 
describe how we (Rajan Gupta, Greg Kilcup and I) do the calculation. For more 
details, see Refs. ^if- We begin by expressing Bk as a ratio 



B 



K 



{K{p = 0)| 37^(l+75)rf sj,{l+j,)d\K{p = 0)) 
0)|s7475f/|0) {0\s^a5d\Kip = 0)) 



imp 



(77) 



Each of the operators in this expression is a shorthand for the lattice operator 
which results from matching with the continuum. Figure |n| shows schematically 
how we did the calculation in Ref. [^. Since then we have changed the method 
slightly, but not essentially |^^. The vertical and horizontal directions represent, 
respectively, space (3-dimensional in practice) and Euclidean time. The expecta- 
tion values indicate the functional integral over configurations weighted, in QQCD, 
by the gauge action. In practice this means an average over some number of con- 
figurations generated with the correct measure. The lines with arrows are quark 
propagators, and the boxes represent the lattice bilinear and quadrilinear operators. 
Finally, the wavy lines at the edge are "wall sources". These create the quark (or 
antiquark) with equal amplitude across the entire timeslice, and thus ensure that 
the mesons which are created have p = 0. Gauge invariance is maintained by first 
fixing the source timeslices to Coulomb gauge. The boundary conditions are peri- 
odic in space, but Dirichlet in time, so that the propagating quarks can bounce off 
the ends of the box, but not propagate through them. 
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Figure 14: Schematic depiction of the method used to calculate Bk- 



We use wall sources because they give us the freedom to insert any operator 
we wish for the "boxes" in the diagram. For example, the operators can involve 
quarks and antiquarks at slightly different lattice sites. This freedom turns out to 
be crucial for staggered fermions, because (as explained below) the operators we 
want to use involve q and q at different positions. More generally, it allows us to 
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Figure 15: Data for Bk (with tree- level matching). The quark mass is such that the 
lattice kaon is slightly heavier than the physical kaon. Sample of 23 configurations. 



test that different discretizations of the continuum operator yield the same results. 
The disadvantage of wall sources is that (in the way we implement them [Q) they 



create not only kaons, but also i^*'s and other excited strange mesons. These give 
contributions which fall off like exp[— (m^* — mft')r], where r is the distance from the 
source. Thus if one is far enough away from both sources only the kaon contributes 
to the matrix element. I give an example of our data in Fig. This shows the 



ratio of Eq. as a function of the timeslice on which the operator resides. The 
operator has been summed over all space, which significantly improves the signal, 
ans is another advantage of the wall sources. The desired signal is independent of 
r, since for all r a particle of mass mx (either a i^o or a Kq) propagates the length 
of the lattice. Edge effects due to excited states and particles bouncing off the 
boundary are apparent, but there is a "plateau" covering a considerable number 
of timeslices from which to extract the signal. We improve the statistics further 
by averaging over this plateau, although, since the results at different times are 
correlated, the improvement is not as significant as one would naively expect. 

To give you a feel for how the analysis proceeds, I show in Fig. |1^ the results 
at four lattice spacings, and for a variety of values of rriK- These are for operators 
matched to the continuum only at tree level. The kaon mass has been converted to 



physical units using the hadron spectrum to set the scale |Q . The dashed vertical 
line shows the value of the physical kaon mass. 

The statistical errors in the results are small, in part due to the use of a ratio 
to calculate Bk- The small errors allow us to clearly observe the following features: 

• Bx has a smooth dependence on m|-, and, apparently, a finite chiral limit. 
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Figure 16: Results for using tree-level matched operators. The lattice spacings 
are roughly 0.2, 0.105, 0.08 and 0.06 fm as one descends the plot. 



• There is no need to extrapolate to get to the physical kaon mass. Quenched 
particles with the kaon mass can be simulated directly on present lattices. Of 
course, this is a cheat, since the results in the figure are mostly for degenerate 
quarks, and one must extrapolate to the non-degenerate case. 

• There is a clear and significant dependence on lattice spacing. 



43 



6.2 Matching continuum and lattice operators 

Lattice operators are chosen so that, at tree level, they have the same matrix 
elements as the continuum operators, for momenta much lower than the cut-off 
[p -C vr/a). With Wilson fermions it is easy to find such operators. A continuum 
bilinear or four-fermion operator is discretized into a lattice operator having exactly 
the same form, with all fermion fields on the same lattice site. The tree-level matrix 
elements are the same as those of the continuum operator for all the momenta that 
are available on the lattice. 

The matrix elements of lattice and continuum operators do not, however, agree 



when one includes loops. Examples of one-loop diagrams are shown in Fig. [17. 
Lattice propagators and vertices differ from their continuum counterparts when 
the momentum in the loop is of 0(1/ a). We have already seen this difference for 
fermion propagators — compare Eq. ^ with the continuum propagator. It is true 
also for gluon propagators {D~^ = 4 sin(/i;^/2)^ versus fc^), and the quark-gluon 
vertex. To do the one-loop matching, one must add parts to the lattice operator, 
proportional to g^, so as to make the matrix elements agree. These additional terms 
are finite, because they come from short distances and the lattice integrals have an 
ultraviolet cut-off. Since the momenta involved are ~ vr/a, the corrections should 
be reliably calculable using perturbation theory, as long as a is small enough. 

Let me mention a subtle point that one tends to forget when doing the matching 
calculations. It is not sufficient that perturbation theory be reliable at the scale vr/a. 
One actually needs the stronger condition that it be reliable down to (0.5 — I) /a. 
This is because, to do the matching, one must, in principle, compare matrix elements 
with external Euclidean momenta satisfying \p\ ^ Aqcd- This is necessary to avoid 
the infra-red region where perturbation theory breaks down. But one also must have 
pa small enough that lattice artifacts in the matrix elements are small. Typically 
this remains true for pa < 0.5 — 1. In practice, at one-loop one uses a gluon mass 
to regulate the infra-red divergences, and sets the external momenta to zero. This 
is adequate, because in the matching calculation, the infra-red contributions are 



identical and cancel. This point is emphasized in Ref. [^|. 

One-loop continuum calculations are straightforward. The corresponding lat- 
tice integrations are, however, a mess, and are evaluated numerically. I can per- 
sonally attest that, soon after beginning such a calculation, one asks oneself the 
questions: Is this matching really necessary? Can't we just use the lattice regu- 
larization throughout? Unfortunately, the answers are yes and no, at least for the 
moment. The point is that the electroweak theory has a chiral representation of 
fermions, so its discretization is problematic. The chirality is not an obstacle to 
discretizing a left-handed operator such as that in (Eq. |T^), because only the 
even-parity part of the operator contributes, which is the average of left-handed 
and right-handed operators. Incidentally, it is possible to directly match from the 
electroweak theory including weak bosons to the lattice effective Hamiltonian, and 
then run down to a low scale on the lattice. 
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Figure 17: Diagrams contributing to one-loop matrix elements needed for matching. 



The result of one-loop matching takes the general form 



lat 



0{g') + 0{a). (78) 



Here z is a set of continuum operators which mix with each other under the con- 
tinuum RG. To define these operators, which in general contain 75, one must pick 
a renormalization scheme as well as a renormalization scale. One of the standard 
schemes is naive dimensional regularization (NDR). The set of lattice operators 
which are required for matching is, in general, larger than the set of continuum 
operators, [jj] Thus the finite matrix Cij is rectangular. The anomalous dimen- 
sion matrix is, however, square. It governs the dependence of the continuum 
operators on /i. 

Let me give a concrete example, relevant to the present subject. In order to 
simplify the allowed Wick contractions consider four-fermion operators composed 
of four distinct flavors 



Oll = i [Vi7/.(l+75)^2 ^37/.(l+75)V'4 + (2^4) 

S = i [V^iV^2 V^3V^4 + (2 ^ 4/ 

= I [V'i75^2 V'375V'4 + (2 ^ 4)] , etc. 
The result of matching for Wilson fermions is 



Qcont 
LL 



1 + 



9' 



167r2 
2 



-41n(^) 

TT 



54.753) 



+ [csS + CpV + ctT + c,V + CaA] + 0{g' 



(79) 
(80) 
(81) 



^2) 



^3) 



"The appearance of extra operators is not peculiar to the lattice. Even in the continuum, when 
one uses dimensional regularization, extra "evanescent" operators are needed at intermediate stages 
of the calculation, associated with the additional — 2e dimensions. 
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The most important feature of this result is the appearance of lattice operators 
having all possible tensor structures. In the continuum, the matrix elements of Oll 
are constrained by chiral symmetry to vanish as mj^ in the chiral limit. This is 
not true for the matrix elements of Os, O-p, etc, which couple to both LH and RH 
quarks. Q This problem is due to the breaking of chiral symmetry by the Wilson 
fermion action — even though the breaking is 0(a) at tree- level, divergent loops 
give factors of 1/a leading to finite contributions at one-loop. The consequence 
is that if matching is not done exactly, to all orders in g"^, matrix elements will 
not have the correct continuum chiral behavior. In addition, there will be lattice 
artifacts, proportional to a, with the wrong chiral behavior. Here we have the 
doubling problem coming back to haunt us. This proves to be a significant obstacle 
in practice, and, because of this, it is preferable to calculate Bk using staggered 
fermions. Probably the only hope for similar accuracy with Wilson fermions is to 
use non-perturbative matching . 



I want to make two general comments about matching. As we have seen, lattice 
calculations must be combined with continuum coefficient functions obtained from 
RG equations. To be consistent, if one uses 1-loop matching, one must use 2-loop 
RG equations. To see this, recall the solution to the 2-loop RG equation for the 
case of an operator which does not mix. Running from a heavy scale rriH down to 
yU, one finds 



c(/i) = c{mH) 





7(0)/2/3o r 









167r2 ^2/3o 2/32 

(84) 

where 7*^"^ and Pn are the (ri + l)-loop contributions to the anomalous dimension 
and /^-function, respectively. The term in the rightmost parenthesis proportional to 
g^ifi) is of the same form as a one-loop matching correction. To be consistent, one 
must include this term in c(/x). Since it contains 7*^^^ and it comes from two-loop 
running. 

My second comment is that the perturbative matching can be done equally well 
in QCD and QQCD. Indeed, at one-loop the two theories give the same results for 
fermionic operators, since fermion loops do not enter. The issue does arise, however, 
when picking the scale at which to evaluate c(yu) . The product of c(/i) calculated in 
QCD, with BxifJ^) evaluated in QQCD, is not independent of /i because /3o, Pi and 
7i receive contributions from quark loops and thus differ in the two theories. One 
must simply guess a value of fi for which it seems reasonable that QQCD will do 
the best job of imitating QCD. 



**The only exception is V + which has the same positive parity part as Oll, and does have 
vanishing matrix elements in the chiral limit. 
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6.3 Staggered fermions and chiral behavior 



I now give a lightning review of the essentials of staggered fermions. For more 
details see the two texts, or my articles explaining in detail how and why one uses 



staggered fermions to calculate matrix elements involving PGB[^ 
We begin with naive fermions 

nfi n 



and perform a change of variables, known as "spin-diagonalization" [^O 
Note that 7„ depends only on mod2(n^). The result is 



S 



N 



E 



^5) 



^6) 



^7) 



The gamma matrices have been replaced by the phases (thus the name of the 
transformation) 

r?/.(n) = (-)^"<."''. (88) 

Note that the four spinor components of x ^ire only coupled by gluon exchanges. 
Staggered fermions result from simply deleting three of the four components of 
X, leaving a one component fermion on each site. This reduces the number of 
continuum fermions from 16 to 4. 

Staggered fermions represent four degenerate flavors in the continuum limit, but 
at finite lattice spacing the SU{4) fiavor symmetry is broken. The only continuous 
symmetry of the action Eq. is that corresponding to fermion number 



X{n) exp[i6'y]x(n) , x(^) ^ x(^) exp[-z6'y] . 



^9) 



In the massless limit, the continuum symmetry enlarges to SU{'i)i x SU{4)fi, while 
the lattice action has an additional symmetry 



XW 



exp [i [ 



^OA]x{n) , x{n) x{n)ew[i{-)"d, 



(90) 



where (— )" = (— )"i+"2+™3+™4_ This turns out to be a fiavor non-singlet axial sym- 
metry in the continuum limit, and we refer to it as f/(l)A- It is very important 
as it guarantees that mass is multiplicatively and not additively renormalized (so 
that setting the lattice parameter m to zero really means massless quarks). It also 
guarantees that Bx has the correct chiral behavior [|5^, ^ ^ . The remaining sym- 
metries are translations, rotations, refiections and charge conjugation — the discrete 
residue of flavor and Poincare symmetries. 

To use staggered fermions, we must learn how to identify the continuum spin 
and flavor transformation properties of flelds constructed from x and x- To do this 
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requires some notation. We denote quark fields in a continuum theory with four 
degenerate fermions using upper case letters, e.g. Qa,a- These have a spinor index 
(here a) and a flavor index (here a), both running from 1 to 4. Quark bilinears in 
this theory are of the general form 

Qa,alfeFQpfi ^ Q{lS®iF)Q. (91) 

75 determines the spin of the bilinear, the flavor. The Dirac matrices are labeled 
using the notation of Eq. ^ except that now is one of the 16 "hypercube" 
vectors having components which are either or 1. It is convenient to label the 
possible flavors in a similar way (rather than using the generators of the group 
S'f/(4)), using the complex conjugate matrices = 7f [H ■ 



With staggered fermions, the discretization of the continuum fleld Q is spread 
over the 16 positions of a 2'^ hypercube[^. Divide the lattice into such hypercubes 
(choosing one of the sixteen alternatives) , labeled by the positions of their corners, 
y = (0,0,0,0), (2,0,0,0), (2,2,0,0), etc. Label positions within hypercubes by a 
hypercube vector A^. Then the 4x4 matrix fleld 

XQ(l/)/3,fe = |Ex(2/ + ^)(7a)/3,6 , (92) 

A 

corresponds to the continuum fleld Q as a ^ 0. In particular, the discretization of 
the continuum bilinears is obtained by substituting xq ^or Q in Eq. In general, 
these involve x and x flelds at different sites on the underlying lattice, and so must 
be made gauge invariant. This is usually done either by joining quark and antiquark 
by an appropriate string of link matrices, or by flxing to a smooth gauge, such as 
Landau gauge. In the new notation the axial symmetry is 

Xq exp [2^a(75 ® ^5)] Xq , (93) 

which shows that it is a flavor non-singlet. 

To write the action in this new notation requires extending it to include deriva- 
tives. There are various choices. To study the corrections due to discretization, 
it is convenient to keep the invariance of the action under translations by a single 
lattice spacing as transparent as possible. A notation with this feature is given in 
Ref. |5^. The action is simply 



ExQ(y)[(7M ® ^)D, + m]xQ{y) , (94) 
y 

which looks just like the continuum action. The complications of the lattice are 
hidden in the deflnition of D^, which mixes up xiv) with the fleld in adjacent 
hypercubes. For this reason the action is S'f/(4) invariant only in the continuum 
limit. The recovery of this invariance can be made more explicit using a form of the 
action written entirely in terms of hypercube flelds , but I do not display this 



result here as it will not be needed in the ensuing discussion. 
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We are now in a position to understand how to find a lattice operator which 
matches onto the continuum operator needed to calculate Bx, i.e. 



cont 



(95) 



(color indices now explicit). Staggered fermions show their dark side in such a 
matching calculation. We choose to introduce a different staggered fermions for 
each continuum quark. This means that the lattice theory does not become QCD 
in the continuum limit, but rather a related theory in which there are four versions of 
each quark. Thus the matching really proceeds in two steps [^. (1) Find operators 
in the extended theory whose matrix elements involve the same contractions as does 
Ocont between a physical K and K. In this step we divide out the effect of the extra 
fermions by hand. (2) Find the lattice operators which match onto the operators 
in the extended theory needed in step (1). 

Step (1) can be done in many ways. Our choice is to use the following operator 
in the extended theory (for a different choice see Ref. IST 



O' 



cont 

A' 
A'' 



A' + A'' + V' + V'' 

Sa{l^.l5 ® i^)Db 5;(7^75 ® i^)D'^ 



16 



(96) 
(97) 
(98) 
(99) 
(100) 



This must be sandwiched between an initial kaon created by the source D' {'-^^®^^)S' , 
and a final K destroyed by D{j^ ® ^^)S. The primes are a device to restrict Wick 
contractions so that the spin and fiavor indices are contracted in two-loops, as 
sketched in Fig. 11^. This is necessary because Fierz transformations are completely 



different in the four-fiavor theory from those in QCD. It is simple to convince oneself 
that the Wick contractions of 0'^^^^ in this particular matrix element are the same 
as those of Ocont between a physical K and K, the factors of 16 canceling the 
contributions from the extra quarks. 




Figure 18: The two spinor loop contraction of the vector operator. 



To carry out this first matching step, we could just as well use any other fiavor 
in place of ^5, as long as we change the fiavor in both the operator and the 
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external sources. Our choice of flavor ,^5 is only important once we have matched 
onto lattice operators. The point is that the U{1)a symmetry has flavor ^5, so that 
the lattice PGB arising from chiral symmetry breaking has this flavor. Only the 
matrix elements involving this pion satisfy the Ward Identities which are needed to 
show that the chiral behavior is the same as for the continuum operator . 



Step (2) must be done perturbatively. At tree level, the result is an operator of 
exactly the same form, written in terms of the lattice flelds xs, Xd, Xs' and xd'- At 
higher orders, flavor symmetry breaking in the lattice action means that matching 
introduces many of the 65,536 operators of the form 

Xsils ® ^f)xd Xs'ils' ® iF>)XD ! (101) 

Fortunately, at one- loop, we need only consider the subset of these operators with 
flavor ^5. The contributions of operators with other flavors are suppressed either by 
powers of a ox g"^. Higher loop calculations would likely be a nightmare, as many 
more operators must be included at intermediate stages. 

The one-loop result for operators made gauge-invariant by flxing to Landau 
gauge is 



57, 55 



0;„„.(NDR.ri = Ol(i + ^[41„^-1i) (102) 

+ [lo.T^L, + 1.7 - lo.ivL - 7.5v,2l . (103) 

lOTT^ L J 

The calculations have also been done for two other types of operator: gauge invariant 
2^ operators with gauge links joining the x aiid x flelds [^, and "smeared" Landau 
gauge operators, which are spread across a 4^ hypercube in such a way that their 
tree- level matrix elements have no terms of 0(a)[Q. The corrections turn out to be 
small, particularly for the 2^ and 4^ Landau gauge operators. For the former, this is 
apparent from the above result. Recalling that (7^/(167r^) 1/80, typical corrections 
are seen to be ~ 10%. In fact, there are cancellations, so that the overall correction 
is only a few percent. The corrections are nevertheless signiflcant, as they are larger 
than the statistical errors. This is illustrated in Fig. ffUl which shows an example 



of the difference between tree- level and one-loop matching for PS | . 

When using the matching equation, one must choose a value of g^, i.e. choose 
both a scheme and a scale. (Note that the scale in need not he the same as 
II.) When making these choices one is, in effect, guessing the higher order terms. 
As discussed above, a physical coupling scheme such as MS leads, in general, to 
more rapid convergence. As for the scale, a reasonable choice would be to use the 



BLM prescriptionj^^ (the simpler scheme of Ref. [TJ] does not apply to a matching 
calculation). Unfortunately, this requires knowledge of the parts of the two- loop 
graphs involving fermion loops, which have not yet been calculated. So one is left 
to make reasonable guesses. We make two choices: /i = vr/a, the lattice cut-off, and 
yU = 2GeV , the scale at which we quote our result. We use the average of the results 
from these two choices as our central value, and half the range as an estimate of 
the systematic error. The resulting error is small, ~ 2%. 
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Figure 19: Tree-level [g"^ = 0) vs. one-loop {g^ = 1.82) matching for Landau gauge 
2^ (unsmeared) and 4"^ (smeared) operators. 



6.4 Chiral perturbation theory 

ChPT allows one to predict the expected behavior of B^- At leading order, the 
operator Ocont (Eq- iH) is represented in the effective theory of PGB's by 

|/^5(Ea^St),,(S9^St),, . (104) 

(For the notation see, for example, the lectures of de Raphael.) / is the pion 
decay constant at leading order (in the normalization that = 93 MeV), while 
B is the leading order value for Bk- ChPT does not give any information on B, 
but does determine the form of the chiral expansion, and the normalization of the 
non-analytic contributions from chiral loops. One finds (with rriu = m(i)||25|| 



Bk = B 



g2 

l-i3 + -)y\ny + by + cye'' + 0iy^) 



(105) 



where y = mj^ / {Att f)"^ ^ 0.2 is the usual chiral expansion parameter, while e = 
{nis — rrid) / (m^ + m^) measures the non- degeneracy of the s and d quarks, b and c 
are unknown constants. The important features of this result are 

• Bx is predicted to have a finite chiral limit. This means (given the definition 
Eq. |1|) that the matrix element of Ocont vanishes as mj^ in the chiral limit. 
This can be traced back to the fact that the operator contains only left- 
handed fields. It is this behavior which is maintained on the lattice when 
using staggered fermions, but not with Wilson fermions. 
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The dependence on e is non-leading, suppressed by ylny.Q Thus there should 
not be a large change in when changing from a theory in which ms{lat) = 
ma{lat) = ms{phys)/2, to one in which nisilat) = ms{phys) and md{lat) ^ 0. 
A rough estimate is given by assuming the chiral log gives the entire effect 

0.03 -0.05. (106) 

The range comes from varying / between and fx, and varying the cut-off 
in the log from rrip to 47r/^. 

To really calculate this ratio one needs to know c, one of the large number of 
higher order constants that appear when one extends the chiral Lagrangian 
to include the electroweak effective Hamiltonian. One might have hoped that 
c could be measured from another weak decay, but this turns out not to be 
the case . 



What happens in the quenched approximation? The result becomes 



l-{3 + e^)y\ny + b^^^y + c^^^ye^ + 6\'^—^\n 



2e VI + e 



(107) 

All except the last term is contained in Ref. [p5|. The term proportional to b is 



a new result [^. The most noteworthy feature of Eq. |107| is that, for degenerate 



quarks, the form is identical to that for QCD. (The term proportional to 5 vanishes 
when e = 0.) In particular, the chiral logarithmic term has the same coefficient 
in both theories. Because of this, I would argue that the QQCD result for Bk is 
likely to be more reliable than that for quantities for which the chiral logs differ 
between QCD and QQCD, e.g. f^, and ra^j f^,. This does not necessarily mean 
that Bk will be the same in QQCD and QCD, because the constants i?(°\ h^^^ and 
c^o) are different from their QCD counterparts. But, in some sort of average sense, 
quantities for which the chiral logs are same in the two theories are likely to be 
more similar. How it is that the chiral logs are the same for Bk is explained in Ref. 

The QQCD and QCD results differ, however, for non-degenerate quarks. Indeed, 
the term multiplying 5 diverges in the limit e ^ 1. This is similar to the divergence 
noted above in m^. Thus there is no sense in using slightly non-degenerate quarks in 
QQCD and trying to extrapolate to e = 1: the 5 contribution is a quenched artifact, 
and the logarithmic term proportional to differs in the two theories. Thus I think 
that the best approach is use QQCD to study Bk for degenerate quarks, and wait 
until we have QCD itself under control before studying the physical kaon. For QCD 
an extrapolation from small e may be possible. 

^^That the dependence is quadratic in the mass difference foUows from the "CPS" symmetry 
introduced in Ref. IpOl. 
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6.5 Errors due to quenching 

The only way to really know the effect of quenching is to calculate Bk in QCD. An 
important step in this direction has been made by two recent calculations which 
include the fermion determinant |]63|, [6^. Both calculate -Bx(e = 0) for theories with 



two degenerate quarks, the masses of which range from down to ms/2, Although 
these theories are a fair distance from QCD — their "pion" clouds are composed of 
particles with the masses of kaons — these results do provide very useful information 
about Bk in QCD. 

These calculations are done on 16'^ lattices at a ~ Ifm, roughly the smallest 
volume and lattice spacing that it is reasonable to use. This illustrates one way in 
which QQCD will remain useful as QCD calculations improve. We can use quenched 
calculations to study the dependence of quantities on the volume, lattice spacing 
and quark masses, and map out the regions in parameter space where the systematic 
errors are small. We then do initial QCD calculations at the edge of the allowable 
parameter space, and use the quenched results to estimate systematic errors. Later 
on we check these errors directly in QCD. 

The results of the two simulations are simple and striking: they agree with 
quenched results at the same lattice spacing, within ~ 5% errors. This agreement 
is both for the slope of Bk vs mj^ and for the intercept. 

What is the significance of this result? First let me sound a note of caution. This 
is a result at only one lattice spacing. As shown in Fig. |T^, Bk has a substantial 
dependence on the lattice spacing in QQCD. The agreement at a ~ 0.1 fm could be 
fortuitous, with the lattice spacing dependence being different in QCD and QQCD. 

Assuming this is not so, let us analyze the result in terms of the chiral expansions 



Eqs. |105| and |107| . The result for two flavors, and degenerate quarks, has exactly 
the same form as for QCD and QQCD [^. Thus what the two simulations show 
is that 

and 6(2) ^^{o)_ ^^Qg) 
Extrapolating linearly in the number of flavors, this suggests that 

B = fi(^) ^ 5^°) and b = b'-^^ ^ 6^°) , (109) 

in which case 

BK{e = 0)^BP{e = 0). (110) 

This result is exactly what the hand-waving of the previous section suggested, but 
it is perhaps surprising how small the difference between QCD and QQCD is. Of 
course, it needs to be checked by adding a third quark. 

What these simulations do not address is the value of c — the coefficient which 
determines the dependence of Bk on — m^. As explained above, the chiral logs 
suggest that the extrapolation to e = 1 will raise Bk by ~ 5%. 
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6.6 Symanzik's improvement program 

The final source of systematic error comes from working at finite lattice spacing. 
If we calculate at a variety of different spacings, but always match to the same 
continuum renormalization point, we expect 

BK{^DR,fi)a = BK{NBR,fi),ont[l + aAi + a^Al + 0{a'')] . (Ill) 

The An are non-perturbative scales characterizing the discretization errors. In fact, 
we will see that the first correction term is absent (Ai = 0) if the calculation is 
done with staggered fermions in the manner described above. This is in contrast to 
Wilson fermions, for which the 0{a) term is present, and diverges like l/mj^ in the 
chiral limit. 

Let me sketch how one arrives at Eq. Ill using perturbation theory [^|. One 



compares matrix elements of continuum and lattice operators between quark states, 
for some set of external momenta, but now keeping terms proportional to powers of 
a. At tree level, powers of a come only from the discretization of the operator. At 
higher orders, one also has to worry about corrections to propagators and vertices. 
For example, the staggered fermion propagator involves 

sin(A;^a) = ak^{l - + ...), (112) 

(no implicit sum on /i). Note that there are only even powers of a; the same is 
true for the vertices, but not, in general for the operators. For Wilson fermions the 
corrections to propagators, vertices, and operators are all, in general, of 0(a). 

At 0{g^), one can show, by power counting, that the 0{a) terms can multiply log- 
arithmically divergent integrals. Thus one gets corrections of the form pag"^ In(am), 
where p is some combination of the external momenta or masses, and m a similar 
combination providing the infra-red cut-off to the integral. Since oc l/ln(a) at 
small a, this term is no smaller than tree-level contributions. Similarly n— loop 
terms proportional to ag'^"' In(am)" are of the same size. The conclusion is that the 
entire 0(a) term in Eq. 111| is of the form aAi, where Ai is a non-perturbative scale. 



A similar result holds for the a^ and higher order terms. 

Symanzik has laid out a systematic program for removing discretization errors 
ip. It has been proved to work for scalar theories, and there is no reason to doubt 
that it works also for gauge theories. I will provide only a thumbnail sketch. There 
are two parts. 

1. Add all possible dimension 5 terms to the action, 

S4 — > 5'4 + a^Q(5f^)S'5i , (113) 



and adjust their coefficients so that some set of physical quantities does not 
have 0(a) corrections. This is to be done order by order in perturbation 

^n=0 ' 



theory, q = ^^"^^=0 '^ind'^"' ■ At tree-level, one adjusts qo so as to cancel the 
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0{a) terms in S4. The strong result of Symanzik is that this removes terms 
of the form a{g'^lna)^ for aU n, at least for physical quantities. What remains 
are terms of size ag"^, which are removed by determining the Cji, etc. A very 
important constraint is that the S^i have the same symmetries as S4. This 
is because one is trying to cancel lattice artifacts due to 5*4, and these are of 
restricted form due to translation, rotations, parity and any other symmetries. 

2. If one wants to "improve" the matrix elements of operators, one must not only 
improve the action, but also improve the operators themselves. For example, 
for the dim-6 operators appearing in B^, the improvement is made by 

Oe c{g')0, + d,ig')Oy, . (114) 

j 

Tree-level improvement removes a(5'^ln(a))" terms, leaving ag"^, etc. The op- 
erators that one has to add are all those with the same symmetries as Oq. 
If this is a local operator it will have fewer symmetries than 5*4, and thus 
typically there will be more possibilities for Ojj than for 5*5^. 



I now sketch the application of this program to staggered fermions 35, 56 . As 



noted above, the staggered fermion Lagrangian can be written in a compact form 

-54 = E x(2/)(7m ® i)D^x{y) + mx{y){i ® i)x{y) ■ (115) 

y 

The advantage of this notation is that the only bilinears which transform as singlets 
under translations are those with = I- This and other symmetries restrict the 
possible dimension 5 operators to be 

X{I®I)D^D^X, x{cr^u(E)I)F^,x and m^xil I)x ■ (116) 

All these three terms are, however, forbidden by the U{1)a symmetry. (To see this 
the mass must be treated as a spur ion field which rotates under the symmetry.) 
There are thus no dimension 5 operators available to improve the action. This in 
turn must mean that there are no 0{a) corrections to on-shell quantities, because, 
if there were, then from the nature of the improvement program, they could be 
removed by the addition of 5*5. This is true to all orders in perturbation theory, 
and I shall assume it also holds non-perturbatively. 

Now I turn to the operators needed for the lattice calculation of Bk- Let Cf^^ 
be a vector composed of all lattice operators which appear when we match to the 
continuum operator we are interested in. As discussed above, this is a long vector. 
Let Cconi be the corresponding vector of continuum operators. The dimension 7 
operators that we must subtract from Of^^ will be all those with the same sym- 
metries. The symmetry group is much smaller now because these are quasi-local 
operators; translations and rotations reduce to those transformations which map 
the 2^ hypercube into itself — the hypercubic group. The transformation properties 
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of staggered fermion bilinears under this group has been studied by Verstegen ^7 . 



The net result is that one ends up with an extremely long vector of dimension 7 
operators, which I call Oj^^. I now assume that Symanzik's results apply, so that 
it is possible to find an improved lattice operator whose matrix elements, in lat- 
tice states, agree with those of the continuum operators in continuum states, up to 
corrections of O(a^) 

{c{g')Ol + d{g')aOl) = {OlJ{l + 0{a')) . (117) 

Here c{g'^),d{g'^) (square and rectangular matrices, respectively) are perturbatively 
calculable coefficients, obtained by removing 0{a) terms in a set of matrix elements. 
The claim is that this improvement then applies to all other matrix elements. 

To demonstrate the desired result I need to relate Oj^^ to its continuum counter- 
part. But dimension 7 lattice operators mix with dimension 6 operators when one 
includes loops — the ultraviolet divergences give contributions oc 1/a which cancel 
the a which comes with the operator .0 Thus to obtain a continuum dimension 7 
operator one needs to subtract out some dimension 6 lattice operators. I assume 
the form 

{d{g')Ol + d{g')aOl) = {aOlJ{l + 0{a)) . (118) 

Notice that this is only good up to corrections of 0(a) and not 0(a^)n. Combining 
Eqs. |117| and |118| one finds 



{OD = -rT^((OL* - dd-'aOlJ){l + Oia')) . (119) 

c — da '-c 

The product of matrices multiplying the matrix element on the r.h.s. is just a 
perturbative matching matrix, part of which has been calculated to one-loop, as 
described above. 

The final part of the argument is messy but straightforward. One simply enu- 
merates the operators Oj^^, using the results of Refs. [^, |5^. An example is 



X5'E(7/.75®e5)^MXD' XsiI®Y.^-)XD. (120) 

The continuum version has the same form, but with Xs ~^ S, etc. The important 
result is that in none of the operators do both bilinears have fiavor ^5. Thus, when 
we insert the operators between states of flavor ^5, their matrix elements vanish 
identically in the continuum, where flavor is a good symmetry. In other words 

{OLh, = ^AiOlJU^ + 0{a')) . (121) 

c — da 

Thus the corrections to all of the lattice dimension 6, AS* = 2 operators in these 
matrix elements are of O(a^). This is true whatever order in perturbation theory one 

ttThere would be mixing with yet lower dimension operators, if any were available, but to make 
a AS' = 2 operator requires at least four-fermion fields. 
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Figure 20: Extrapolating quenched to the continuum hmit. 

does the matching. Truncating the matching calculation at one-loop just gives rise 
to errors of 0{g^), and not of 0{ag'^). This result greatly simplifies the extrapolation 
to the continuum limit, since one knows not to include linear terms in a. 1 do not 
know whether there are terms of O(a^), though I suspect that they can be ruled 
out by similar considerations. 

With Wilson fermions one actually has to do the calculations, rather than just 
discuss formal manipulations. Fortunately, there are many fewer operators because 
there are no additional flavors. The action, together with various interesting op- 
erators, have been improved at tree level, and the one-loop matching coefficients 
have been calculated. This reduces the discretization errors to 0{g'^a). The higher 
order terms in perturbation theory coming from tadpole diagrams have also been 
calculated, which appears to remove the bulk of the 0{ag'^"') terms. 

6.7 Status of results for 

In Fig. show how our results extrapolate to the continuum limit [^. Based on 
the discussion of the previous subsection, we assume quadratic extrapolation — the 
data itself is not good enough to differentiate quadratic from linear dependence on 
a. (The results at a ~ 0.2 fm are too far from the continuum limit to be usefully 
included in the fit.) We find that the scale characterizing discretization errors is 
quite large, A2 ~ 1 GeV. For this graph, we have fixed the lattice spacing using 
rrip, and taken g^ in the MS scheme at /i = vr/a. The results for the two different 
types of operator disagree slightly at the largest lattice spacing, a ~ 0.1 fm, but are 
consistent when extrapolated. This is as expected: two different operators will, in 
general, have different discretization errors. To quote a result we average that for 
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the two operators and use half the difference as an estimate of the systematic error. 

To estimate the other systematic error in our extrapolated result, we repeat 
the calculation using f^, to determine a, and repeat it using matching coefficients 
evaluated with the coupling scaled down to n = 2 GeV. There are also some small 
systematic errors, having to do with our numerical method, which I will not discuss. 
The final result is 

Si^(NDR,2GeV) = 0.606 ± 0.020 (statistical) (122) 

± 0.014 (choice of g^) (123) 

± 0.009 (choice of lattice spacing) (124) 

± 0.004 (choice of operator) . (125) 

In the past I have combined the systematic errors in quadrature, but this is not really 
justified. What is important and striking is that the statistical errors are small, and 
the systematic errors smaller still. Let me also mention that the Japanese staggered 
group have also calculated Bk with 2^ Landau gauge and gauge-invariant operators. 



at (3 = Q and 6.3 and find results consistent with those shown above [£3[. So the 
numerics appear under control. 

This result is for degenerate quarks in the quenched approximation. The error in 
Bk due to these approximations is not known, but, as discussed above, is likely to 
be small, conservatively ~ 10%. 



6.8 Other matrix elements from QQCD 

Many other weak matrix elements are being calculated using lattice methods, al- 
though the calculations are not as advanced as that for Bx- I have no time to 
discuss these in any detail. We have been attempting for years to calculate the ma- 
trix elements leading to the CP-conserving and CP-violating kaon decay amplitudes 

{K\Oi-s\Tr7r) . (126) 

Oi are the operators appearing in the effective electroweak Hamiltonian. Of par- 
ticular interest are Oqs, which give the largest contribution to e'. Since there are 
two particles in the final state, these calculations are done indirectly. ChPT is used 
to relate the matrix element to the simpler, and doable, {K\0\7t) and {K\O\0). To 
raise the level of these calculations to anything approaching that of Bk will require 
either considerably more more statistics and/or an improved method. 

There are by now extensive results for various decay constants (/^, fx, • • • ) 
and for the semileptonic form factors K — > neu, D (vr, K)ev. These have been 
calculated successfully using Wilson fermions, since chiral symmetry does not give 
important constraints. The statistical precision has reached the 10% level, and the 
perturbative corrections are understood. What remains to be done is the extrapo- 
lation to the continuum limit, and for some estimate of the quenching errors to be 
made. 
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7 Heavy mesons on the lattice 



B mesons will be the major focus of experimental studies of CP violation in the 
next decade. Some measurable quantities are directly related to elements of the 
CKM matrix, without hadronic uncertainties. Others are not, and here the lattice 
can contribute. Examples include 

• fs and Bb, which allow an extraction of Vtd and Vts from measurements of 

— Bd and Bs — Bg mixing. 

• The form factors for B {D^D*)ev (which in the heavy quark limit are 
related to the Isgur- Wise function) . These will allow an extraction of Vbc- 

• The form factors for B (7r,p)ei/, which will allow an extraction of Vhu- 

• The matrix element for B K*'~f, allowing extraction of Vts- 

• The quark mass mj,. This is a fundamental parameter of QCD, and an input 
into (or prediction of) models of physics beyond the standard model. 

The hope is that not all the results for CKM elements will agree, thus giving us an 
indication of new physics. 

My discussion in this section will focus mainly on theoretical issues. Numerical 
results have not yet settled down — for example there are still significant discrepan- 
cies between results for Jb in the limit that the b is static. For an excellent review 



of the status of numerical results see Ref. [|33 |. 

An important point to keep in mind is that the properties of hadrons containing 
one heavy quark are determined by non-perturbative, long distance physics, just 
as for light quark hadrons. In other words, it is no easier to calculate J'b than 
/tt. Similarly, the problems introduced by quenching have the same character as 
for light quark hadrons — 5-mesons have PGB clouds, which are absent in QQCD, 
and are replaced by an rj' cloud. Quenched chiral perturbation theory for heavy- 



light mesons has been developed by Booth [^, and by Zhang and I [^. Typical 
estimates of the effect of quenching on decays constants are 10 — 20%. 

There are, however, two simplifications. The first is that the constraints of chiral 
symmetry, central for kaons, are much less important for 5-mesons. Thus Wilson 
fermions are the discretization of choice for the light quark. The complications 
associated with staggered fermions are probably not worth the effort, though not 
all agree with this assessment ! 



The second and most important simplification concerns the dynamics of the b- 
quark. As discussed in the lectures of Grinstein, when ttiq — > oo, the heavy quark 
travels at constant (Minkowski) velocity, and acts as a spinless color-triplet source. 
Furthermore, quark and anti-quark dynamics are decoupled. These observations are 
formalized in heavy quark effective theory (HQET), a methodology which allows 
one to systematically expand quantities in powers of l/mq. Standard examples of 
the results of HQET are 
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TUB* = ruB + 0{m^^] 



The large mass behavior of the decay constant of a heavy-hght meson P in 
which the 6-quark is replaced by a fictitious quark with a different mass 



a{mp) 


2//3o 


_a{mB)_ 


= 0OO 
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B 



+ ... 



m. 



(127) 



The four form-factors involved in B 
function, the Isgur-Wise function ^(f 



(D, D*)eh' are all related to a single 



What the lattice can do is to calculate the unknown functions in these expressions, 
and to check the importance of the non-leading corrections. 

To calculate 0oo5 or ^{v-v'), or other quantities which are present in the uiq — > oo 
limit, one must discretize the HQET. There are some subleties associated with the 
continuation to Euclidean space, but these appear to be understood [^. Things 
are simplest when the heavy quark is at rest. As mentioned in the discussion of the 
confining potential, the propagator for such a static quark is proportional to the 
product of time directed U matrices 



G{n, t; m, 0) = 5^ 



^m,0),4 f^(m,l),4 • • • 



(m,t-l),4 



(128) 



When combined with the usual light antiquark propagator, this makes a heavy 
meson propagator from which one can extract 0oo- The fact that the heavy quark 
does not move increases the statistical noise — a light quark propagator wanders all 
over the lattice and thus averages out some of the fluctuations in the gauge fields. 
To reduce the noise one has to work hard at producing good sources. 

In order to estimate the size of the corrections to the heavy quark limit, one 
must put "dynamical" heavy quarks (i.e. quarks moving at a variable velocity) 
directly on the lattice. There is a potential obstacle to doing so: at present lattice 
spacings the 6-quark Compton wavelength exceeds the lattice spacing, m;,a ~ 2 — 3. 
6-quark propagation thus includes large lattice artifacts. Two approaches have been 
adopted to remove these artifacts. 

The first is to work at small enough quark masses that lattice artifacts are small, 
and then extrapolate up to the b quark. Crudely speaking this requires m^a < 1, 
implying ruq < 3.5GeV for present lattices. Using an improved action probably 
increases the maximum mass that can be used. This method has been used for 
Jb, B ^ K*'y and the Isgur-Wise function. I show in Fig. |2l| results for 0p from 
Ref. [|7^. The point at 1/Mp = comes from a calculation with a static heavy 



quark, while the other points come from dynamical (Wilson) heavy quarks. These 
are results from f3 = 6.3, corresponding to 1/a ~ 3.2 GeV. The fit to the form of 
Eq. |127| is good, and one can read off fpy/m^B from the intercept with the vertical 
line marked "5" . 

This is not quite the whole story. The two points with the smallest non-zero 
values of 1/Mp, have mqa ~ Mpa > 1, and are likely to be afflicted by substantial 
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lattice artifacts. Indeed, an approximate way of removing these artifacts has been 
used to correct all except the static point, and the corrections are substantial (as 
large as ~ 50%) for the points at small 1/Mp. I explain the source of this correction 
below. The bulk of it is reasonable, but there remains a systematic uncertainty in 
these large mass points. Fortunately, this uncertainty does not have much impact 
on fs, because the curve is pinned down from both sides of the B mass. Even if 
were to discard the two points with the smallest values of 1/Mp, the result for Jb 
would be similar. The present status has been summarized by Soni |7^] , who quotes 
fs = 173 ± 40 MeV (in the normalization where /,r = 135 MeV). 

The alternative approach is to use yet another effective theory to describe the 
6-quark, namely non-relativistic QCD (NRQCD) [^. If uiq » Aqcd, then in its 
couplings to low momentum gluons, the quark is non-relativistic. The couplings to 
high-momentum gluons must be treated relativistically, but these can be accounted 
for using perturbation theory. The NRQCD Lagrangian is 

+ ci{g^)-^^^a ■ Bij + c2(/)-^V^t<T ■ {D x E - E x D)^ (129) 



where is a two-component field. The coefficients Cj are obtained by perturbative 
matching with QCD. At leading order they are all unity. 

The crucial point is that, when one discretizes NRQCD, the errors are no longer 
determined by rriQa, but rather by pa, where p ~ Aqcd is a typical momentum. 
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Thus, in heavy-light systems, discretization errors for heavy and Ught quarks are 
comparable. In the the bb system the typical momentum is somewhat higher, p ~ 
Mv ~ O.IM ^ O.SGeV, but still satisfies pa < 1 for present lattices. 

NRQCD teaches us something about heavy Wilson quarks. The point is that, 
as long as uiq » Aqcd, the quarks are non-relativisitic, and must be describable 
by a non-relativistic effective Lagrangian. This will have the same form as Eq. p.29| 
(since this contains all terms), but the effect of discretization errors will be to change 
the coefficients Cj substantially from those in NRQCD |[73|| . The coefficients can be 
determined using perturbation theory. By suitably changing the normalization of 
the fields one can get the leading term in £nrqcd with the correct normalization, 
while by changing the definition of tuq one can obtain the second term up to pertur- 
bative corrections. After these corrections, results for heavy Wilson quarks should 
only have errors of 0(Aqcd/'"^q) (from the incorrectly normalized third term in 
^nrqcd), and not of O^mqa). These are the corrections that have been applied to 



points in Fig. 21 



Nevertheless, the action is not that of NRQCD. Kronfeld and Mackenzie have 
suggested studying heavy quarks using a modified Wilson action |]7^. Their scheme 
interpolates between the standard Wilson fermion action for light quarks and NRQCD 
for heavy quarks. For my purposes here, however, there is no important distinction 
between this and the NRQCD approach, and I will focus on the latter, for which 
more results are available. 

NRQCD is a non-renormalizable theory, because >Cnrqcd contains operators of 
dim-5 and higher. It must be formulated with a UV cut-off, which is here provided 
by the lattice, A ^ 1/a. As mentioned above in the discussion of B^, loop dia- 
grams mix the higher dimension operators with those of lower dimension. Thus, for 
example, there are contributions to wave-function renormalization proportional to 
F{g'^)A/mQ = F {g"^) / {mga) . These are calculable order by order in perturbation 
theory. Since in practice one can only work to finite order, it is clear that one cannot 
take a too small, for then the uncalculated corrections will become large. This is 
not a practical problem at present lattice spacings. 

This does bring up an important issue that, in my view, remains to be fully 
resolved. Even if one could calculate F{g'^) to all orders in perturbation theory, there 
might be non-perturbative terms proportional to exp[— n/(2/3o(?^)] oc (Aqcdo)"" (/?o is 
the first coefficient in the QCD /9-function) . If there are such terms with n = 1, they 
give non-perturbative contributions to wave- function renormalization oc Kqcd/M, 
which are not calculable. This in turn means that there is an ambiguity in the 
A/Mp contribution to (pp (Eq. |127| ) of size oc Aqcd/^p- But A ~ Aqcd, so the 
ambiguity is of the same size as the quantity we wish to extract. This argument 
was first given in Ref. |j7^ . 

There is no dispute over whether such terms can exist (examples are infra-red 
renormalon ambiguities) , but what is controversial is their size [|7B| . These are non- 
perturbative effects at short-distances, and the lore is that these are small. This is 
true of the contributions of short distance instantons, for example. How can this 
dispute be resolved? One can think of these terms as non-perturbative contributions 
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to the matching between QCD and NRQCD. Thus one should compare the results of 
perturbative matching to those of non-perturbative matching — the latter requiring 
either physical matrix elements, or quark and gluon states in a particular gauge. 
Such a program has been initiated in Ref. HSl. An alternative comparison is 



obtained by calculating physical quantities in different ways. If the answers agree, 
the non-perturbative ambiguities are likely small. The example of rrih is discussed 
below. I think it is important to study this issue further. 

The aim of the NRQCD program is to directly simulate 6-quarks. The first step 
in this program is to make sure that one obtains a good description of the bb system. 



The ordering of the terms in Eq. |129| is according to decreasing importance in the 
bb system [^. Terms in the first line are of 0{mQv'^)^ where v is the velocity of the 
heavy quark, and give the spin-averaged splittings (e.g. the IP-IS splitting discussed 
a long way above). The second and third lines contain terms of 0{mQv'^), the former 
being the leading contribution to the hyperfine splittings, the latter correcting the 
spin-averaged splittings. Since v"^ ~ 0.1, the first line gives fine structure to 10%, 
the next line hyperfine structure to 10%, and the third improves the accuracy of 
fine structure to 1%. 



All these terms have now been included in the simulations [p9| , [77[] . Propagators 



can be calculated in a single pass, because the time derivative can be discretized 
as a forward difference. Thus the calculations are much faster than for light quark 
propagators. The status of the calculation of the spin-averaged spectrum is shown 
in Fig. NRQCD refers to Ref. fT^], UK-NRQCD to Ref. 0. Notice that 
simulations with two (moderately light) fiavors of dynamical fermions have now 
been done. The spectrum is little changed with Uf = 2 from that in QQCD; both 
are in good agreement with the experimental data. The lattice spacing has been 
adjusted to fit the IS" — IP and 2S — 2P splittings (from which as can be determined 
as discussed above). Reasonable agreement for spin splittings is also obtained. 

What about m^? Spin-averaged splittings are quite insensitive to nib; a reason- 
able value has been used in the simulations. An accurate determination of rrib can 
be made in two ways. First, one can use mya = 2mf,a + -Enrc^ ~ 2EQa, i.e. adjust 
mfe until twice its mass plus the binding energy (determined from the simulation) 
equals the physical T mass. The measured binding energy must be corrected for the 
self energy of the quark [EQa). This correction can be calculated in perturbation 
theory, E^a = Cig^ + 0{g'^). It is an example of a quantity which might contain the 
previously discussed uncalculable non-perturbative terms oc Aqcd^- These would 
lead to an uncertainty in the extracted value of rrib of size Aqcd- 

The second method is to measure the kinetic energy of the bb states, e.g. 

aEr{p)=aEr,^K+ ^''f +.... (130) 

One adjusts the bare lattice mass m°a until aMr,kin agrees with the physical mass 
in lattice units, aM^. One then uses perturbation theory to match the lattice 
mass to the pole mass, nib = {1 + c'g^ + 0{g^))ni^. In this method an unknown 
non-perturbative term is truly a small correction. 
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Figure 22: Spin-averaged T spectrum from (a) NRQCD (uf = 0): filled circles; (b) 
NRQCD {uf = 2): open circles; and (c) UK-NRQCD (nj = 0): boxes. Experimental 
results are the dashed horizontal lines. 

Both methods lead to consistent results, with errors of ~ 200MeV. Thus a 
difference between the two methods of approximately Aqcdj which is what would 
be expected from non-perturbative terms, is not ruled out. More accurate data is 
needed to resolve the above mentioned dispute. Combining the results from the two 
methods, the final quoted values for the pole mass are 



The next stage is to apply these methods directly to B-mesons. This is just 
beginning. It is my hope that the issue of non-perturbative ambiguities can be 
resolved, and the lattice can, in due course, give results for a number of transition 
amplitudes including I /Mb corrections. 



I hope to have convinced you that lattice calculations are now reliable enough to 
make significant contributions to phenomenology. As time progresses, the impor- 
tance of these calculations will increase. To give an idea of what might happen I 
have played the following game. The lattice results for and fsVBB constrain, 
respectively, Im{V^j) and \Vtd\^ (ignoring small charm quark contributions). Recall- 
ing that Vtd = AX^{1 — p — 17]), we can convert these into constraints on p and rj. I 
want to imagine how these constraints might look, say, five years hence. 

It is my guess that by that time we will have made enough progress with simu- 
lating QCD that the error in will be roughly the same as the present error in 



mbijif = 0) 
mb{nf = 2) 



4.94±0.15GeV 
5.0±0.2GeV. 



(131) 
(132) 
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Figure 23: Possible future constraints on p and rj. 

the quenched result, Eq. |125| . For purposes of illustration, I will also take the same 
central value. I think that the errors in and Bb will be reduced substantially, 
and I assume fsy/B^ = 200 ± 10 MeV, instead of the present 200 ± 40 MeV. Let me 
further assume that the experimental errors in Vcb and Vub will drop significantly. I 
take \ Vcb[= 0.038, iKb/Kfel = 0.080±0.004, and xa = 0.67±0.02 (this is the measure 
of B — B mixing). Finally, I assume rrit = 172 GeV. The resulting constraints on p 
and 7] are shown in Fig. The parameters are supposed to lie between the three 
pairs of curves (the hyperbolae are from K — K mixing, the large circles from B — B 
mixing, and the small circles from |K(fe/V^fe|). The reduction in errors has moved 
the members of each pair much closer than at present, and there is no solution for 
p and 7]] It will be most interesting to see how things develop in reality. 
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